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ABSTRACT 


The  problem  of  reflection  and  refraction  of  weak  hydromagnetic 
disturbances  at  surfaces  of  discontlniiity  in  the  medium  is  studied. 
It  is  shown  that  the  introduction  of  an  appropriate  set  of  direction 
angles  for  the  unperturbed  magnetic  field  on  each  side  of  the 
discontinuity  makes  possible  (l)  a  simplification  of  V.  C.  A. 
Ferraro's  laws  relating  to  the  reflection  and  refraction  of  Alfv^n 
waves  and  (2)  a  simple  graphical  method  for  determining  the  real 
angles  and  speeds  at  which  the  reflected  and  refracted  waves  emerge. 
The  relationship  of  the  graphical  procedure  to  Huyghens  '  constiniction 
is  discussed.   Application  of  the  theory  is  made  to  several 
two-dimensional  steady-flow  problems  and  to  the  problems  of  reflection 
and  refraction  at  a  rigid  infinite  conductor  and  at  a  contact  discontinuity. 
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I.   IMTRODUCTICM 

In.  the  first  paper  of  this,  series  the  propagation  of  weak  hydro- 
magnetic  dlstiirbances  In  a  medium  free  of  "boundaries  and  discontinuities 
was  discussed.   The  medium  was  assumed  to  be  a  classical  hydromagnetic 
fluid  -  i.e.^  infinitely  conducting^  neutral^  perfect  and  compressible. 
The  object  of  the  present  work  is  to  extend  this  analysis  by  discussing 
the  reflection  and  refraction  of  weak  hydromagnetic  disturbances  at 
boundaries  and  surfaces  of  discontinuity  in  the  medium. 

The  problem  to  be  studied  may  be  formulated  as  follows:  Given  (l) 
a  surface  of  discontinuity /J  separating  different  states  of  the  undisturbed 
mediimi,  (2)  the  initial  position  at  an  Incident  wave  front^  and  (5)  the 
strength  of  the  disturbance  at  this  wave  front;  find  the  position  and 
strengt"  of  the  scattered  wave  fronts  at  all  later  times.   The  problem 
falls  naturally  into  two  parts:  (a)  the  derivation  and  analysis  of  the 
hydromagnetic  analog  of  Snell's  laws  and  (b)  the  calculation  of  the 
strengths  of  the  scattered  waves  in  terms  of  the  strength  of  the 
incident  wave. 

The  first  problem  is  complicated  by  the  fact  that  the  disturbance 
speeds  appearing  in  the  laws  of  reflection  and  refraction  depend  upon  the 
direction  of  propagation.   The  particular  form  of  this  dependence  makes 
it  impossible  to  obtain  explicit  formulas  for  the  angles  of  reflection  and 
refraction  in  terms  of  the  angle  of  Incidence,  imless  only  Alfven  waves  are  in- 
volved or  \mless  special  orientations  of  the  magnetic  field  are  assumed.  We  shall 


J.  Bazer  and  0.  Fleischman  ,    "Propagation  of  weak  hydromagnetic  discontinuities". 
Research  Report  No.  MH-10,  Institute  of  Mathematical  Sciences,  New  York 
University,  New  York  (1959)-   See  also  J.  Fluid  Mech. ,  2,  366  (1959). 


present  here  a  new  analysis  of  the  laws  of  reflection  and  refraction. 
This  analysis  is  based  on  the  introduction  of  an  appropriate  set  of 
angles  for  specifying  the  orientation  of  the  magnetic  field  with  respect 
to  the  plane  of  incidence,  specifically,  the  angles  \lr„  and  9„  of  Figure 
1.   The  use  of  these  angles  has  the  basic  advantage  of  laying  bare  the 
symmetries  Inherent  in  the  laws  of  reflection  and  refraction.   This 
enables  one  to  arrive  at  a  relatively  complete  qualitative  picture  of 
the  dependence  of  the  angles  of  reflection  and  refraction  on  the  orientation 

of  the  magnetic  .field  and  on  other  parameters  of  the  problem.   As  a 

2 
by-product,  one  obtains  a  simplification  of  V.  C.  A.  Ferrarc's  laws 

determining  the  angles  of  reflected  and  refracted  Alfven  waves.   More 

important  still, the  use  of  these  angles  makes  possible  a  simple  graphical 

3 
procedure  for  determining  the  angles  and  associated  disturbance  speeds 

of  the  scattered  waves  Induced  by  a  given  incident  wave.   Briefly 

described,  the  procedure  involves  the  intersection  of  a  straight  line 

in  the  plane  of  incidence  with  the  curves  obtained  by  intersecting  the 

surface  of  wave  normals  with  the  plane  of  Incidence.   Furthermore,  when 

applied  to  linearized  steady-flow  problems  [_see,  e.g..  Section  Vj  the 

method  yields  for  each  Mach  number  a  clear  picture  of  the  dependence  of 

the  hyperbolic,  elliptic  and  " hyper- llptlc    flow  regions  on  the 

remaining  parameters  of  the  problem.   The  intimate  relation  of  the  method 

to  Huyghens '  classical  construction  will  be  discussed  later.   Here, 

it  should  be  mentioned  that,  among  other  advantages,  the  method 

presented  seems  better  suited  to  numerical  work  than  graphical  procedures 


^.  C.  A.  Ferraro,  Astrophys.  J.  U$_   (1954). 
^More  precisely,  real  angles. 

This  term  was  introduced  by  J.  E.  McCune  and  E.  L.  Resler,  Jr.  in  their 
study  of  magnetozerodynamic  flows  past  thin  bodies,  J.  Aero/Space  Sci., 
27.  h93   (i960). 
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based  on  Huyghens '  construction)  it  gives  more  accurate  results  and 
furnishes  more  direct  information. 

The  problem  of  obtaining  explicit  formulas  for  the  strengths  of 
the  scattered  waves  is  complicated  by  the  fact  that  each  type  incident 
wave  will,  in  general^  give  rise  to  three  reflected  and,  when  transmission 
is  possible^  three  refracted  waves.   In  calculating  these  coefficients 
one  is  led  to  evaluate  6x6   or  at  best  5x3  determinants.   Since  the 
elements  of  these  determinants  are  themselves  complicated  expressions^ 
the  evaluation  becomes  excessively  laborious  even  in  the  5^3  case. 
Moreover^  the  final  results  are  difficult  to  assess  unless  one  is  willing 
to  assume  special  directions  for  the  incident  wave  [e.g.,  normal  incidence^ 
and  for  the  prevailing  magnetic  field.   Accordingly^  we  have  selected, 
from  the  very  outset,  problems  for  study  which  admit  of  simple  explicit 
solution  by  virtue  either  of  a  special  choice  in  the  direction  of 
incidence  or  of  the  magnetic  field.   In  all  cases,  the  magnetic  field 
is  assumed  to  lie  in  the  plane  of  incidence  to  take  advantage  of  the 
decoupling  of  the  Alfven  waves  from  the  slow  and  fast  waves.   In  the 
first  problem,  we  discuss  the  reflection  of  slow  fast  and  Alfven  waves  by 
a  rigid  infinite  conductor  to  illustrate  among  other  things  the  phenomena 
of  multiple  reflection  and  the  production  of  a  pair  of  genuine  hydromagnetic 
waves  [slow  and  fast^  by  a  single  soun.d  wave  incident  along  the  magnetic 
field.   In  the  second  problem,  we  study  the  scattering  of  nonnally 
incident  slow  or  fast  waves  by  a  contact  discontinuity.   In  this  case  it 
is  shown  (l)  that  the  boimdary  conditions  can  be  met  by  a  pair  of  slow  or 
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fast  scattered  waves  [according  as  the  incident  wave  Is  slow  or  fast]  . 
and  that  the  "distiirbance  energy"  is  conserved  in  the  scattering  process. 
Finally,  we  apply  our  resiilts  to  a  series  of  problems  relating  to 
two-dimensional  steady  flow.  We  show,  for  example,  how  to  calciilate  the 
lift  due  to  pressure  forces  on  a  thin  wing  moving  steadily,  with  small 
angle  of  attack,  along  the  magnetic  field. 

The  reflection  and  refraction  of  hydromagnetic  waves  at  various 

types  of  interfaces  has  been  discussed  by  several  investigators  among 

2  5 

whom  are  V.  C.  A.  Ferraro,   and  P.  H..  Roberts.   We  shall  discuss  the 

points  of  contact  of  these  works  with  our  own  later  on,  in  the  appropriate 

context.   Here,  it  is  enough  to  say  that  both  authors  assume  incompressi- 

bility  and  employ  planar  disturbances  which  vary  sinusoidally  in  space  and 

time. 

For  the  sake  of  simplicity,  it  will  be  assumed  here  and  in  the 
subsequent  discussion  that  the  state  of  the  unperturbed  medium  on  each 
side  of  the  surface  of  discontinuity  oU  is  everywhere  constant  and  that 
J^  and  the  incident  and  reflected  wave  fronts  are  planar.   It  should  be 
stressed,  however,  that  the  entire  discussion  applies  locally,  in  the 
neighborhood  of  the  common  curve  of  intersection  of  txT^  and  the  incident 
wave,  when  Ju    or  the  incident  wave  is  curved. 

II.  THE  LAWS  OF  REFLECTION  MD  REFRACTION 
A.  Preliminary  Remarks 

Figures  1  and  2  depict  the  basic  geometrical  setup.  (S'^represents 
the  surface  of  discontinuity  separating  the  upper  half-space  ^  from  the 


^  P.  H.  Roberts,  Astrophys.  J.,  121,  720  (1955). 
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♦  A^V 


FIGURE    2 


lower  h.alf-space^ .  The  vector  H  is  the  prevailing  magnetic  field 
at  the  origin  0.   The  local  plane  of  Incidence  at  the  origin  0, 
determined  by  the  normal  on  the  incoming  wave  front  and  the  imit  normal 
N  on  (^  at  0^  Is  assumed  to  coincide  with  the  (x^z)-plane  and  is  denoted 
by  (^ .   As  shown,  ^   Is  directed  along  the  positive  z-axis.  The  vector 
H  is  the  projection  of  H  on  the  plane  of  incidence.   The  quantities 
)k„.   0  <  \lf„  <  jt/2  and  9^^.  -jt  <  Q„  <  jc,  are  the  angles  between  H  and  H 
and  H  and  H ,   respectively.  We  shall  have  occasion  to  refer  to  the  angles 
3  and  7  of  Figure  1  later.   It  has  been  implicitly  assumed  in  Figure  1 
that  jH  is  continuous  across  o©*  .   This  need  not  be  the  casej  e^'^may^for 
example, be  a  shock  front.   In  general,  if  Q  is  any  quantity  in  2^ ,  we 
shall  employ  Q  to  denote  the  corresponding  quantity  in  "^  whenever  a 
difference  exists. 

Figure  2(a)  shows  the  wave  fronts  in  the  plane  of  incidence. 

W  denotes  the  incoming  wave  front,  propagating  toward  ^a^  f rom  below. 

W.  is  a  fast,  AlfvSn  or  slow  wave  front  according  as  1  =  1,  2,  or  J.  ^j-,: 

W   and  W   denote,  in  the  order  stated,  the  reflected  fast,  Alfven  and 

slow  wave  fronts  and  W.,  ,  W^^  and  W.^  the  refracted  slow,  Alfven  and  fast 

14   15      lo     

wave  fronts  -  all  produced  by  the  Incident  wave  of  type  (i).  Figure  2(b) 
shows  the  disposition  of  the  noimals  n.  and  n. .  on  these  wave  fronts; 
to  keep  the  diagram  uncluttered  only  three  normals  have  actually  been  drawn. 
The  angle  0. ,.-rt/2  <  9  <  rt/2,  is  the  angle  between  N  and  the  incident 
wave  normal  n^.  For  j  =  1,  2,  3,  Q^y    [©^j  <  ^   when  9^  >  Oj  9^^  >  -it  when 
9  <  ol  denotes  the  angle  between  N  and  the  j-th  reflected  wave  -  normal 
n  .5  for  j  =  *+,  5;  6,  9   ,  -rt/2  <  6.  .  <  3t/2;  denotes  the  angle  between 
N  and  the  j-th  transmitted  wave. 

The  n^   and  n^^'s  are  chosen  to  point  along  the  direction  of  propagation 
of  the  wave  front. 
'^'Angles  In  the  plane  of  incidence  measured  from  the  positive  direction  of 
N  Ire  regarded  as  positive  if  their  sense  Is  clockvd.se,  otherwise  negative. 
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In  the  plane  of  incidence  (see  Figure  2(b))  the  incident  and  scattered 
waves  intersect  in  a  common  point  0  which,  it  will  he  supposed,  moves 
with  the  V  =  Vx  at  the  instant  depicted.   Here,  x  denotes  the  uxdt 
vector  along  the  positive  x-axis.  V  will  hereafter  be  referred  to  as  the 
trace  velocity.   Let  U  >  0  be  the  speed  at  0  in  the  normal  direction  n. 
of  the  incident  wave  and  let  U. .  >  0,  j  =  1,  2,  . . . ,  6  denote  the  local 
normal  speeds  in  the  directions  n^^  .  of  the  scattered  waves.   Then, as  in 
geometrical  optics  and  acoustics,  it  is  easily  shown  [See   Appendix  ij 
that  (l)  the  n. .  lie  in  the  plane  of  incidence  and  (2)  for  any  i 


%=U,   , 


(2.1) 


V  .  n,  .  =  U,  .  , 

^      -ij    ij 


j=l,2,...,6    (2.2) 


In  cont.-ast  to  geometrical  optics,  however,  the  U  . 's  are  direction-dependent. 

Q 

The  U  's  are,  in  fact,  given  by 


U.  =  u  ♦  n^  +  c. 
1   "^p   ~ii  —  1 


i  =  1,  2,  5 


(2.5) 


where^ 


2        2      2 
(l+r)  -  kr   cos  ^if^   cos  (O,  -  9„) 
Jl        In 


,   (2.i^) 


~  =  cos  tjj  I  COS  (Qg  -0jj) 


(2.5) 


cf.  the  paper  in  reference  1,  equations  (5«5)  -  (5*8) 
9 
The  existence  of  three  disturbance  speeds  was  discovered  by  N.  Herlofson, 

Nature,  I65,  1020  (195O)  and  independently  by  H.  C.  van  de  Hulst  in  a 
paper  which  appeared  in  Problems  of  Cosmical  Aerodynamics  (Central  Air 
Documents  Office,  Dayton,  Ohio  I95I) ,  p.  k6. 
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2        2      2 
(l+r)  -  kr   cos  \(r„  cos  (O-.-  9^) 
n.  ^    n 


^  1 

2 


1 
2 


.  (2.6) 


In  these  expressions,  u  is  the  component  of  the  unperturbed  fluid 

***p 

2—1  1  /p 
velocity  in  the  plane  at  incidence:  b  =  {[lE   p   )  '   is  the  Alfven 

disturbance  speed  in  the  direction  of  the  magnetic  field  Hj  p  is  the 

density  in  ??  j  r  =  a  /b  is  the  squared  ratio  of  the  sound  speed  a 

and  the  Alfven  speed  bj  and  c  ,   c  and  c  denote,  in  the  order  stated, 

the  slow,  Alfven  and  fast  disturbance  speeds.   It  follows  directly  from 

equations  (2.5)  -  (2.7)  that  (l)  c  /b  is  an  increasing  function  of  \(f  , 

while  c^ /b  and  c  /b  are  decreasing  functions  of  \|r  ,  r,  0  and  0., 

j  =  1,  2,  5  being  assumed  fixed;   (2)  for  fixed  \(r  and  r  and  all  9 

and  e. 


-^  <  min  (l,r)  <  max  (l,r)  <  ^   ; 


(2.7) 


(5)  for  fixed  \|/  and  r  and  all  0  =  0  =  0 


Ci  <  Cg  <  c^ 


(2.8) 


The  corresponding  expressions  for  the  U_,  .  and  c.  .  are 

ij      ij 


U^  .  =  u  .  -  n.  .  +  c.  .  , 
Ij   "v-pj   ~ij  -  ij 


j  =  1,  . .. ,  6, 


(2.9) 


and 


^ij 


(l  +  r.)   -  li-r  .cos  i,^.cos  (0.  .  -  0tt-) 


,   J=i^6, 


(2.10) 


c.  . 

-^     =     COS  ^iTgj    |cos(9^j   -   Q^.)\    ,  j   =  2,5.  (2.11) 


ii     =         __Ii  _  ^   [fl+r.)^  -   hv.    cosV.cos^fQ.  .-  9„.)r 


J 


^     (l+r.)^  -   kv.    GOs\„.cos^(9,  .-  9^.)  r 


j   =   5,^.  (2.12) 


Here,  u    .   is  defined  by 

.wpj 


u    . 


r^      >      ^    =   1.2,3, 


lu^     .     J  =  ^.5.6, 


(2.15) 


and  the  quantities  r.,  b  . ,  \|r   and  9   are  defined  in  an  analogous 

manner.   In  equation  (2.15)  the  primed  quantities  refer  to  the  region 

cK.    3XiA   the  unprimed  quantities  to  the  region  ^T-   From  equation  (2.9) 

it  follows  that  in  a  given  direction  n. .,  two  scattered  waves  of  type 

j  are  possible  whenever  lu  .  *  n . . I  exceeds  c. .. 

'-PJ   -*ij  '  ij 

Combining  equations  (2.5)  and  (2.9)  with  (2.1)  and  (2.2),  we 
find  that  for  each  i,  i=l,2,5 

(V  -  u  ,)  •  n,   =  +c.     ,  (2.li4-) 

^«~'   —pi    --i    —  1     -^ 

(V  -  u  .)  •  n, .  =  +c, .     ,   j  =  1,...6,  (2.15) 

or  equivalently  that 

V  |cos(9^-a)  I   =   c^    ,  (2.16) 

.   Y  |cos(9^^-a)|  =   c^^    ,   J  =  1,2,5,  (2.1?) 
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V  |cos(©^^-  a')|   =  c^j    .    J  =  4;5.6  .  (2.18). 

Here,  V=  Iv-u  I  ,     V'=  (V-u'l  and  a  and  a'  are  the  angles  that  the 

relative  velocities  V-u  and  V-u'.  respectively,  make  with  the 
s\irface  normal  N  at  0.   From  equation  (2.17)  >   we  readily  conclude 
that  an  angle  of  reflection  9  .  is  real  only  if  the  component  of 
V-u  along  the  direction  of  the  normal  n.  .  exceeds,  in  absolute  value, 
the  disturbance  speed  along_^  ..   The  same  statement  holds  for  the 
refracted  waves., u'  and  a'  replacing  u  and  a,  respectively. 

Equations  (2.l6)  -  (2.l8)  can  be  reexpressed  as 


|cos(0   -a)  I      |cos(0  -  a)  I 

-T^ =  ^ =  -^   .   J  =  1.2,3,    (2.19) 

ij  i         V 


|cos(9   -  a')      |cos(e  -  a')     1 


c.  .  c. 

IJ  1 


.   J  =  i+,5,6,    (2.20) 


which  have  essentially  the  same  form  as  Snell's  laws  for  the  scattering 
of  acoustic  waves !  Introducing 

p   =-^  >0   ,    p    =  ^  >0  (2.21) 

1       G^  IJ       C^^ 

into  these  equations.,  we  find  that 

p      |cos(0      -a)|      =     p    |cos(9    -a)|      =     -3-        ,      j   =   1,2,3,      (2.22) 

p.  .|cos(9      -  a')  I   =     Pj_|cos(9   -  a')  I    =     -^        ,     J  =  ^,5,6.      (2.23) 


Since  the  (p  .,9. .)  may  be  regarded  as  the  polar  coordinates  of  points 
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in a  (p_,0) -plane  ^  It  follows  directly  from  these  equations  that 
(p^  -^6^  •);  J  =  1^2,3^  all  lie  on  the  line  I   through  (p^^©^)  whose 
equation  is 

i:   p  I  cos  (0- a)  I   =^  =  p.  I  cos  (0.-a)|         (2.21+) 

V 

while  (p  .,<^..),   j  =  ^,5;.6,  alllie  on  the  line  f  '  through  (p.;9  ) 
whose  equation  is 

£':   p  jcos  (Q-a')|  =  "^  =  P  J  cos  (0  -a')j  .     (2.25) 

Y' 

Suppose  now  we  plot  p  .  =  l/c  .(0  .) ,  j  =  1^2^3  against  0  .  on  the  same 
polar  (p^0) -diagram.   The  result  is  evidently  a  family  of  three  curves, 
each  cuxve  consisting  possibly  of  more  than  one  branch.   It  follows 
immedic.,tely  from  equations  (2.21)  and  (2.24)  that  the  points  of 
intersection  of  the  line  i  with  the  curves  of  this  family  satisfy  equation 
(2.22)  and  hence  furnish  the  angles  of  reflection  and  the  reciprocal 
normal  speeds.   Similar  remarks  apply  to  the  refracted  waves. 

B.  Reflection  and  Refraction  at  a  Hydromagnetic  Contact  Discontinuity 
1.   The  Graphical  Method 

This  last  remark  is  the  basis  of  our  graphical  method.   To  gain 
some  additional  insight  into  the  method,  it  is  best  to  discuss  it  in 
relation  to  a  particular  problem.  We  choose  for  this  purpose  one  of  the 
simplest  and  perhaps  most  interesting  problems,  namely,  reflection  and 
refraction  at  a  hydromagnetic  contact  discontinuity.    In  this  case,  it 


Many  of  the  results  for  this  special  case  may  easHj"  be  extended  to  cover 

the  more  general  situation  envisaged  above.   To  see  this,  one  need  only 

rename  the  angles  0.-  a,  0_,  .-a  and  0^^  in  equation  (2.19)    0.,  0.  .  and 
^     i   ^  ij        H     ^       ^   -'z     i   ij 

0„  -  a  respectively  and  note  that  the  resulting  equation  is  identical  with 
n 

the  special  case  of  (2.I9)  considered  in  this  section.   Similar  remarks 

apply  to  equation  (2.20). 
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will  "be  recalled,  only  the  density  and  entropy  may  be  discontinuous 
across  the  surface  oO'  unless  the  normal  component  of  the  magnetic  field 
K   vanishes.  When  H^  vanishes,  a  discontinuity  in  the  fluid  velocity 
u  is  permissible;  however,  for  the  sake  of  simplicity,  it  will  be 
supposed  that  u  is  again  continuous  across  au.      Prom  the  Galilean 
invariance  of  the  basic  jump  relations,  it  follows  that  u^^may,  without 
loss  of  generality,  be  chosen  to  vanish  everywhere.   Hereafter,  it  will 
be  assumed  that  this  is  the  case. 

Let  us  temporarily  focus  our  attention  on  the  problem  of  reflection. 
Here,  it  is  a  question  of  finding  the  intersection  of  the  line  i  =  i(|v|/b) 

pblsinOl   =  — -^=—     =     p.blslnO.I  (2.26) 

|v|/b      ^ 

with  the  family  of  curves  obtained  by  plotting 

(p.,b)   =  i ,  (2.27) 


^Ij 


(c,./lo) 


against  9. .,  j  =  1,2,3  on  the  same  polar  diagram.   The  Alfven  speed  is 

-'-J 

introduced  in  equations  (2.26)  and  (2.27)  to  reduce  the  number  of  curve 
parameters.   These  parameters,  which  appear  in  the  expressions  for  c .  . /h 
(see  equations  (2.9)  -  (2.12)),  are  r,  \j;^  and  0^^.   In  Figures  3(a)  -  (c)  , 
we  show  the  polar  plots  of  equation  (2.27)  for  three  values  of  rj  namely 
r  =  1/2,  (Figure  3(a)),  r  =  1  (Figure  3(b))  and  r  =  l/2  (Figure  3(c)). 
In  all  cases  0„  is  taken  to  be  1+5°  and  \|f^  =  0  .    Figure  3(c)  also 


In  Figure  3j  H  therefore  is  Identical  with  H. 

■=■  '^■'   »Vvp  /WV./ 
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contains  curves  associated  with  the  parameter  values  r  =  2.,   6„  =  ^5° 
and  \|f  =  ij-5°}  we  shall  have  occasion  to  refer  to  these  later.   As  in 
Figure  2.,   all  angles  are  measured  from  the  vertical  N-axls  and  are 
considered  positive  or  negative  according  as  they  are  measured  clockwise 
or  counterclockwise  from  W. 


What 


ever  the  value  of  r,  the  oval-like  curves  (labelled  f  ) 

correspond  to  the  fast  p^  branch j  the  pair  of  solid  lines  (labelled  A  )^ 

perpendicular  to  the  H  -axis,  correspond  to  the  Alfven  p   branch  and 

the  hyperbola- like  curves  (labelled  s  )  correspond  to  the  p  -slow 

branch.   The  pair  of  interrupted  lines,  perpendicular  to  the  H  -axis 

are  the  loci  of  points  satisfying  the  equation. 

1 
pb  cos  \lfjj|cos  (0-0^)1   =   [(l+r)/r]^   •        (2.28) 

These  lines  serve  as  asymptotes  of  the  s  -curves  whether  \(f  =  0  or  othervise. 

The  f  ,   A  and  s  -curves  are  all  symmetric  with  respect  to  the  origin 

and  with  respect  to  the  H  -axis.   It  follows  directly  from  equation  (2.27) 

and  (2.10) -(2.12)  that  when  ilr^r  =  0  the  Alfven  branch  intersects  the 

n 

slow  branch,  the  fast  branch  or  both  according  as  r  exceeds  unity,  is 
less  than  unity  or  equal  to  unity.   All  branches  have  zero  curvature 
on  the  H  -axis  when  r  4  !•  When  r  =  1,  both  the  slow  and  fast  branches 
are  singular  on  the  H  -axis.   They  have  points  in  common  at  pb  =  1, 

0=  0  and  pb=  1,  0=  ©„+  n  and  at  each  of  these  points  are  tangent  to  a 

H  H 

pair  of  lines  that  make  an  angle  of  X  =  +tan  (1/2)  with  the  H  axis. 

[See  Figure  5(^)3  .   By  rotating  the  curves  of  Figures  5(a) ,  3(b)  or  5(c) 

about  the  H  =  H-axis  we  arrive  at  a  family  of  surfaces  in  p-space  known 
P 


-  13  - 

collectively  as  the  surface  of  vave  normals. 

The  curves  of  Figures  5(a)  ;  "('b)  and  (c)  are  typical  of  those 

in  the  parameter  classes  r  >  1^  r  =  1  and  r  <  1^  respectively.   Keeping 

r  and  9  fixed^  and  increasing  ij;^,  we  find  from  equations  (2. 10) -(2. 12) 

and  (2.27)^  for  each  9_,  that  the  points  on  the  f  -curves  move  toward 

the  origin  and  the  points  on  the  A  and  s  -curves  move  away  from  the 

origin.   The  general  appearance  of  these  curves^  however^  remains  the 

same  as  in  the  case  \|/  =  0  [see^  e.g.^  the  curves  of  Figure  3(c)  associated 

with  ^  =  lj.5  J  except  on  and  in  the  neighborhood  of  the  H  -axis, 
•ti  p 

Here,  the  double  or  triple  points  -  i.e.,  the  points  where  the  Alfven 

wave  intersects  the  slow  branch  (r  >  l) ,  the  fast  branch  (r  =  l)  or 

both  (r  =  1)  -  split  up  into  two  and  three  separate  points  respectively. 

Moreover,  in  the  case  r  =  1,  the  curvatures  of  the  slow  and  fast  branches 

on  the  H  axis  are  no  longer  singular.   Finally,  we  note  that  the  sole 

effect  of  varying  0  while  keeping  r  and  \|;  fixed  is  to  rotate  the  f , 

-tl  H 

12 
A  and  s  -curves  rigidly  with  respect  to  the  fixed  axes  along  N  and  x^. 

The  effect  of  varying  the  angular  parameters  \j/  and  9  may  be 

rL       n 

visualized  as  changing  the  orientation  of  the  surface  of  wave  normals 
with  respect  to  the  plane  of  incidence. 

In  Figure  h{^&)  ,   we  have  superimposed  the  line  I   of  equation  (2.26) 

on  the  curves  of  Figure  5(a).   The  symbol  M  introduced  in  the  figure 

is  defined  by  11        ^ 

\     -       b  -pb|sin9,|  ^2-29) 

12 

Rather  than  draw  a  new  set  of  s,  A  and  f-curves  for  each  new  9„,  it  is 

H 

evidently  simpler  in  numerical  work  to  deal  with  one  fixed  set  of  curves 

and  to  rotate  the  N-axis  with  respect  to  this  fixed  set  in  such  a  manner  that 

H  makes  an  angle  9„  with  N. 
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and  is  the  reciprocal  of  the  distance  of  i  from  the  N-axis.  M^  will 
hereafter  "be  referred  to  as  the  Mach  number  to  suggest  the  analogy 
wlth  corresponding  gas  dynamical  quantity.   This  terminology  is  clearly 
appropriate  j  for  an  observer  moving  with  the  trace  velocity  V^  will  see 
a  stationary  configuration  of  wave  fronts  into  which  the  fluid  moves 
with  Mach  number  Ivl/b. 

By  varying  the  angle  of  incidence  and  coionting  the  number  of 
intersections  of  i  with  the  lower  portions  of  the  f    ,   A   ,    s  -curves, 
one  arrives  at  the  following  conclusions:   (l)  For  sufficiently  small 
angles  of  incidence  (and  correspondingly  large  Mach  numbers)  there  are 
three  (real)  angles  of  reflection  ^  0_ ,  0.^  and  0.^.   (2)  To  each 
fast  incident  wave  with  sufficiently  small  |0  |  there  correspond  three 

such  angles  of  reflection.   However,  when  0  is  approximately  70°  ttie 

1^ 
fast  reflected  wave  is  critically  reflected.    (J)  When  the  incident 

wave  is  of  the  slow  or  Alfven  type  then  at  positive  and  negative  angles 

of  incidence  which  are  sufficiently  large  in  absolute  value,  the  fast 

wave  is  critically  reflected.   (U)  At  sufficiently  large  positive  angles 

of  incidence,  Alfve'n  or  slow  waves  can  produce  at  most  two  real 

reflected  waves  -  one  Alfven  and  one  slow  wave.   The  same  is  true  of 

negative  angles  of  incidence  which  are  moderate  in  absolute  value. 

(5)  However,  at  negative  angles  of  Incidence  that  become  successively 


■\hether  or  not  the  associated  modes  of  propagation  are  excited  or  not 
depends  upon  the  bo\:indary  conditions  and  on  the  orientation  of  the 
magnetic  field. 

A  wave  of  type  j  due  to  an  incident  wave  of  type  i  is  said  to  be 
critically  reflected  when  0. .  just  becomes  equal  to  +90  or  -90  . 
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larger  in  absolute  value ^  first  the  Alfven  and  then  the  slow  wave  is 
critically  reflected. 

Figure  4(b)  is  a  plot  of  0   ,  0   and  0   against  the  fast  angle 

of  incidence  for  r  =  0.5^  0„  =  k^°   and  trr  =  0.   It  was  obtained  by 

n  n 

graphical  means  from  an  enlargement  of  Figure  k{a) .      Statements  (l)-(5) 

refer  specifically  to  Figures  i4-(a).   However^  with  obvious  changes^ 

they  apply  for  all  values  of  v,  ^^   and  0^^  the  angles  0„  =  +90°^  +l80° 

n       n  n   —     — 

being  excepted.   Statement  (5)  is  no  longer  valid  in  the  exceptional  cases, 
but  (l)-(l4-)  still  apply  with  minor  changes.   We  leave  to  the  reader  the 
analysis  of  the  special  cases  [see  Figures  10(b)  and  llj . 

To  determine  the  angles  of  refraction  we  must  find  the  Intersection 
of  the  line  £ ' 

1       11 


pb  '  sin0  =  p .  b  '  sin0 


1 


i'°(|v|A.)°»b'°77\' 


1/2    ^  /  ,  /  ^l/2  /o  ^n^ 

T\  '     =  p-  =  (p7pJ  '  (2.30) 

with  the  upper  parts  at  the  curves  obtained  by  plotting 

against  0.  .  on  a  polar  diagram.   Note  that  here  (cf .  equations  (2.25)  ^Jid 
(2.25))  we  have  referred  the  disturbance  speeds  to  the  Alfven  speed  b'  in 
the  region  OJ    and  the  quantity  R  '  =(p'/p)   M^  is  the  Mach  nianber  appro- 
priate to  the  region"^.   Observe,  in  addition,  that  equation  (2.50)  is  of 
the  same  form  as  equation  (2.27)  except  that  in  (2.3)  the  ratios  c. ./^' 
depend  upon  r'  not  r.   But,  If  as  we  shall  here  and  hereafter  assume, 
both  7[  and  ^T^  consist  of  polytropic  gases,  then  we  can  conclude  that 


15 

The  properties  listed  in  the  next  subsection  will  be  found  useful  for 

this  purpose. 
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FIGURE  4(b) 
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^2 

7P 

-     7P'        a'2 
2             2 

(2.32) 


since  H  and  the  pressiire  p  are  continuous  across  a  hydromagnetic 

contact  discontinuity.   It  follows  that  c  , /b  '  is  the  same  function  of 

9.,  as  c. ^ /b  is  of  9. , ,  that  c,^/b'  is  the  same  function  of  9.^  as 
i4     il'         il'       ly  15 

c  /b  is  of  9^^  etc.   Thus^  on  plotting  p  .b  '  versus  9  ^^  j  =  h ,'^ ,(i , 
for  r  =  1/2,  1);  =  0  say^  we  obtain  curves  [see  Figure  5]  which  are 
identical  with  those  in  Figure  i|(a)  .   However^  for  a  given  angle  of 
incidence  the  distance  (M,  ')    of  the  line  £'  from  the  N-axis  differs 
from  that  of  £  by  the  factor  r;  '   =  (p'/p)    . 

Referring  now  to  Figure  5;  we  observe  that  for  sufficiently 
small  angles  of  incidence  or  sufficiently  large  Mach  numbers  the  line 
f  =  i'(Nl  ')  will  intersect  the  fast  branch  (labelled  f ,  ) ^  the  Alfven 
branch  (labelled  A  )  and  the  slow  branch (labelled  sA    exactly  once. 
If  the  incident  wave  is  fast,  then  whatever  the  angle  of  incidence, 
there  will  always  be  a  fast  refracted  wave  unless  p  exceeds  p '  in  which 
case  there  will  exist  a  critical  amgle  of  incidence  at  which  the  fast  wave 
is  totally  reflected.    If  the  incident  wave  is  an  Alfven  or  a  slow 
wave  and  the  angle  of  incidence  is  positive^ then  fast,  Alfven  and 
slow-waves  will  be  totally  reflected  at  successively  larger  angles 
of  incidence.   For  successively  smaller  negative  angles  of  incidence 
we  find  (l)  the  fast  angle  of  refraction  is  -90  degrees,  (2)  i'  intersects 
the  Alfven  branch  and  slow  branches  each  exactly  once  (3)  i'  intersects 
the  Alfven  brajich  twice  and  the  slow  branch  once  and  {\)    £'  intersects 


Since  the  angles  of  refraction  are  necessarily  acute  angles,,  only  those 

intersections  are  counted  that  lie  in  the  region  N  >  0. 

17 

At  this  angle  the  refracted  wave  front  is  normal  to  the  surface. 
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both  the  slow  and  the  AlfVen  tranches  twice.  With  obvious  modifications^ 
these  featvires  of  Figure  5  apply  for  all  r,  ^\f     <   3t/2  and  ©„  4  1  p  ■>  ^ 
and  n.      In  the  exceptional  cases  aU.  of  the  above  discussion  again 
applies  with  minor  changes  apart  from  statements  (3)  and  (h)   which  are 
no  longer  tinae. 

2.  Some  General  Properties  of  the  Angles  of  Reflection  and 
Refraction 

The  angles  of  reflection  and  refraction  enjoy  several 

general  properties  which  are  consequences  of  the  particular  way  in 

which  the  (c.  ./^<j)  '^  depend  upon  9.  .  and  the  parameters  0.  ,  0  ^  \1; 

and  r.   We  wish  to  present  some  of  these  now. '  For  this  purpose  it  is 

convenient  to  write  the  laws  of  reflection  and  refraction  in  the  form 


sin  0. .  sin  0. 
ij  _  i 

c_,  ./h  c./'b 

ij  i' 


sin  0. .     sin  0. 

ij  ^        1 

C./b'     Tl^/^(c./b) 


J  =  1.2,5  (2.35) 


,  J  =  ^.5.6,        (2.3i|) 


1/2    b     f    ,  /  \  1/2 
r\  '     =  -rj    =   (p'/p) 


and  to  employ 


to  express  the  functional  relationship  between  0. .,  0  and  the  other 

-1  Q 

parameters  of  the  problem.    The  following  properties  are  then  more  or 


T75 

It  is  understood  here  that  T)  =  1  when  j  =  1,2  or  5- 
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% 
less  direct  consequences  of  these. equations  and  the  formulas  (2.9)-(2.1l) 

.19 


for  the  disturbance  speeds: 
®ij  &h]  =  @ij  ^^J  ' 


j  l^hJ  -  ^-ij  l^^hJ  ^  ^1 


^,,  Cr-^   =®,,  M    ^  =  :^; 


P. 


j  =  l,3,i^;6,  ^ 


1^ 


(Dij  Kj  ej   =   -®  ..  [9.3  -9j  ,  P^    (2.36) 


rt  -9,  ,  if  9.  >  0 


ii  L-©iiJ   =  -  1   ,  ^  \    ^        J  =  ^'2,5, 
^    ^        '  -(:t  +  9.),  if  9^  <  Oj 


^6 


®ji  [^®ij>  ^'^^  =  ®i  ^    J  =  ^^5,6,  P^ 

\\^,\  <  \\^\  <  \^,e\   '  ^9 

In  these  relations^  parameters  not  shown  are  considered  fixed  and 
unless  specifically  indicated  otheivise  the  subscripts  i  and  j  assume 
the  values  1^2^3  and  l;2j...^6,  respectively.   The  angles  9. .  appearing 
in  equations  P^  and  Pq  are  what  we  shall  hereafter  refer  to  as  the 
acute  angles  of  reflection.   They  are  defined  as  follows: 


"ij 


:  -  9.  .,  if  9.  .  >  0        i  =  1;2,5, 

(2.37) 
•(9^^  +  rt),  if  9.^  <  0      j  =  1,2,5. 
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Recall  that  the  first  index  i  refers  to  the  incident  wave  and  the 

second  index  j  to  the  reflected  or  refracted  waves. 
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Equation  P  expresses  the  fact  that  the  angles  of  reflection  and 

refraction  are  Independent  of  the  direction  of  the  component  of  H_^ 

which  is  perpendicular  to  the  plane  of  Incidence.   Equation  P  states 

that  these  angles  are  also  independent  of  the  sense  of  H  .   Together, 

P  and  P  imply  that  the  reflected  and  transmitted  angles  are 

independent  of  the  sense  of  II_.   Equation  P  states  that  the  slow  and 

fast  angles  of  reflection  and  refraction  are  Invariant  under  the 

transformation  of  r->-r   whenever  the  incident  wave  is  slow  or  fast. 

From  P,  we  conclude  that  if  a  wave  incident  at  the  angle  9.  produces 

scattered  waves  at  the  angles  9. . .  then  a  wave  incident  at  the  angle 

ij 

-9.  will  produce  scattered  waves  at  the  angles  -9. .  only  if  the  direction 
1  .   ij 

of  H  is  reflected  in  the  N-axls  specifically,  only  if  9,,  is  replaced 
"~p  "^  H 

^y  -©XT'   If  ©TT  is  not  so  replaced,  then  according  to  P^.  the  scattered 

n        n  P 

waves  emerge  at  the  new  angles  9f .  =  -  (S) .  . 1-9. ;9„1  ,  which  differ,  in 

general,  from  -9. ..   Only  when  9  is  0  ,  +90  and  l80  do  we  find  that 

•)(■ 
9. .  =  -9. ..   In  other  words,  negative  angles  of  incidence  are  not 

physically  equivalent  to  positive  angles  of  incidence  except  for 

special  orientations  of  the  magnetic  field.   The  next  two  properties  P^ 

and  P  verify  the  reversibility  of  the  scattering  process.   Thus,  for 

example,  P  may  be  Interpreted  as  follows:  If  9. .  is  an  angle  of 

refraction  associated  with  the  angle  of  incidence  9.  ,  then  a  wave  in  /(_, 

incident  at  the  angle  9.  .  will  produce  a  transmitted  wave  in  ('J^  of  type 

1  emerging  with  the  angle  9..   In  the  two  relations  Po  and  P  ,  it  is 

assumed  that  9.  .  's  and  9, . 's  are  real:  these  relations  affirm  a  fact 
ij        ij 

to  be  expected  on  physical  grounds,  namely,  that  the  fast  reflected 

and  refracted  waves  emerge  first,  followed  by  the  reflected  and  refracted 
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Alfv6n  waves  which  are  in  turn  followed  by  the  reflected  and  refracted 
slow  waves  (see  Figure  2(a)). 

3,  Numerical  Considerations  -  Anomalous  Reflection  and  Refraction 
For  sufficiently  large  angles  of  incidence  it  will  he  found 
that^  regardless  of  scale ^  the  intersections  of  the  lines  i(R  )  and 
i'(M  ')   [see  equations  (2.25)  and  (2.29)]  with  the  slow  and  Alfven 
curves  fa3JL  outside  the  range  of  the  graph.   For  such  angles  it  is 
necessary  to  supplement  the  graphical  analysis  with  suitable  approximate 
formulas.  Fortunately,  the  Alfven  curves  are  linear  and  the  slow 
branches  are  almost  linear  in  that  they  approach  the  straight  lines 
of  equation  (2.28)  asymptotically.   These  facts  may  be  employed  to 
obtain  analytical  expressions  for  the  angles  of  reflection  and  refraction 
associated  with  large  angles  of  incidence.   Consider,  for  example, 
the  problem  of  determining  the  angles  at  which  the  reflected  slow  and 
Alfven  waves  emerge.   In  this  case,  the  intersection  of  £   with  the  A 
and  s^-curves  lead  to  the  following  expressions: 

sin  0.„  sin  9. 

i2  1 


cos  tjj  |cos  (©j_2-Qh^I     (Cj_A) 


} 


sin  0., 
i5 


^  \  1/2    sin  9. 


cos  tjj  Icos  (9^^-9g)|  V    /       (c^/b) 

The  first  of  these  equations  is  simply  equation  (2.35)  with  c  /b 
replaced  by  the  explicit  expression  given  in  the  right  member  of  equation 
(2.11)'  this  equation  is  exact.   The  second  equation  is,  on  the  other 
hand,  an  approximate  equation:  it  is  derived  by  replacing,  for  each 
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0, ^,  the  quantity  p   in  the  equation  c   =  p"   by  the  approximating 
value  of  p   on  the  asymptote  [see   equation  (2.28)]  and  then  substi- 
tuting the  resulting  expression  for  c . .,   into  equation  (2.53)«   Such 

an  approximating  value  always  exists  unless  0  =  0^  in  which  case  £ 

20 
is  parallel  to  the  asymptote.    From  the  above  equations  and  the 

trigonometric  s-um  formula  for  the  cosine  it  follows  that 

sec  tjj  sec  0^^.  (c^/b) 
ctn  0^2  =  + tan  0^  , 

sin  0. 

1 

—  1  1  /p 
(l+r"  )  '  sec  f     sec  0  (c./b) 

ctn  0^  ^  =  + tan  0„  . 

^  sm  0 


In  these  equations^  the  "+"  sign  is  intended  when  the  intersection  of  i 
with  the  A- line  or  asymptotic  line  occurs  in  the  half -plane  into  which 
H  is  directed;  otherwise  the  negative  sign  must  be  chosen. 

It  is  worthwhile  noting  the  possibility  of  a  kind  of  reflection 
which  does  not  arise  in  conventional  gas  under  similar  circumstances. 

For  the  sake  of  concreteness^  suppose  that  the  incident  wave  is 
slow  and  the  relevant  Figure  for  the  reflection  process  is  Figure  i4-(a) . 
When  0,,  the  angle  of  incidence^  is  positive  and  sufficiently  small 
the  disposition  of  the  wave  fronts  is  as  shown  in  Figure  2(a).  As  0 

if. 

is  increased  to  "critical  angle",  0^  say  [defined  as  the  angle  at 
which  0   =  90°3 ,   the  fast  wave  front  approaches  the  normal  direction 
to  the  surface  Jo .   What  we  should  like  to  stress  here  is  the  following: 


20 


It  should  be  mentioned  that  it  is  not  difficult  to  obtain  explicit 

formulas  for  0. .  when  0„  =  0. 
ij       H 
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By  increasing  the  angle  9  slightly  past  the  critical  angle  9  ,  it 

is  possible  to  extend  the  manifold  of  fast  wave  solutions  continuously 

past  9   =  90  .   Physically  this  means  that  the  fast  wave  fronts 

associated  with  such  values  of  9-,  emerge  in  the  same  quadrant  as  the 

5 

slow  incident  wave  front  as^  for  example^  the  dashed  line  does  in 
Figure  2(a) .   The  reader  can  convince  himself  from  an  inspection  of 
Figure  4(a)  that  for  negative  9^  one  can  pass  over^  by  continuously 
decreasing  9  .  to  the  case  where  both  the  Alfvln  and  slow  reflected 
waves  emerge  in  the  same  quadrant  as  the  incident  wave.   If  one  were 
to  construct  these  wave  fronts  by  means  of  Huyghens '  method  (see 
Section  B-5) ,   one  would  find  that  the  rays  which  guide  the  "energy 
in  the  mode"  [see  Section  II-B-5  and  III-A-lJ  all  have  components 
directed  along  -N,^  as  do  the  rays  associated  with  ordinary  reflected 
waves.   Incidentally,  in  Figure  \,   the  outward  normals  to  the  s,  A 
and  f-curves,  drawn  at  their  points  of  intersection  with  i   give  the 
directions  of  the  rays  [cf .  reference  l] .   This  "radiation  condition" 
together  with  the  continiiity  considerations  support  the  inclusion 
within  the  manifold  of  reflected  waves  those  which  emerge  in  the 
same  quadrant  as  the  incident  wave. 

The  existence  of  the  above  type  of  anomalous  reflection  process 
[and,  in  fact,  of  others  not  mentioned,  for  the  reason  that  they  are 
less  physicalljr  admissible]  makes  it  necessary  to  define  the  terms 
reflected,  refracted  and  even  incident  wave  more  precisely.   Hereafter, 
we  shall  admit  to  the  class  of  reflected  jjrefractedj  waves  only  those 
which  (l)  may  be  arrived  at  by  continuously  increasing  the  angle  of 
incidence  from  zero,  and  (2)  are  associated  with  ray  vectors  having 


-  25  - 

components  that  are  directed  along  -N  [along  n] .   Finally,  a  wave 
which  propagates  toward  Jj    from  below  -  i.e.,  a  wave  whose  normal 
n  has  a  component  directed  along  N  -  will  be  admitted  as  an  incident 
wave  only  if  the  associated  ray  vector  has  a  component  which  is 
directed  along  N.   It  should  be  mentioned  that  in  Figure  5.  only 
those  admissible  angles  of  reflection  which  are  associated  with 
normally  [as  opposed  to  anomolouslyj  reflected  waves  are  shown. 

Another  interesting  kind  of  anomolous  reflection  (refraction) 
process  is  the  phenomenon  of  conical  reflection  (refraction)  [See 
Reference  ij .   Conical  reflection  (refraction)  can  occur  only  when 
r  =  1  and  H^ is  in  the  plane  of  incidence  and. then  only  at  suitable 
angles  of  incidence.   In  Figure  3(a),  for  example,  conical  reflection 
can  occur  only  when  the  line  I   is  in  the  position  shown  -  that  is, 
only  when  £  passes  through  the  triple  point  where  the  s,  A  and  f 
cuives  intersect.   Imagine  that  the  disturbance  on  the  incident  wave 
front  is  concentrated  in  a  small  disc  of  negligible  diameter.   Then, 
employing  considerations  similar  to  those  given  in  Reference  1  [p.  57^J 
we  can  infer  that  this  "point"  disturbance  will  emerge  as  a  ring 
disturbance  with  an  illximinated  point  in  its  center.   The  disc 
determined  by  the  illuminated  ring  and  its  center  point  is  perpendicular 


to  H  and  moves  parallel  to  H  with  the  speed  b  =  ViaH  p   ;  its  radius 
t  seconds  after  the  incident  wave  has  hit  the  interface  can  be  shown  to 
be  ^t. 

k.   Explicit  Formulas  for  the  Alfven  Waves 

We  assume  here  that  an  AlTven  wave  is  incident  on  the  contact 
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discontinuity^  and  we  derive  explicit  formulas  for  the  angles  at 
whlcli  the  reflected  and  refracted  Alfven  waves  emerge.  We  begin  by 
combining  equations  (2.35)  :>  (2. 5^1-)  and  equations  (2.5)  j>  (2.11)  and 
find,  after  cancelling  the  common  factor  cos  i|/„,  that 


sin  9^2  sin  9^ 


022=03(922-9^)   a^c^os{Q^-Q^)    ' 


sin  92i^  sin  92    .   "  2 


(2.58) 


o^^zos{Q^^-Q^)        o^zos[Q^-Q^) 


(2.59) 


In  these  equations  the  a's  are  +1  or  -1  according  as  the  neighboring 
cosine  factors  are  positive  or  negative.   Employing  the  standard 
trigonometric  formulas  for  the  cosine  of  a  difference  of  angles,  we 
find  easily  that 

ctn  9„^  =  -^^   (ctn  9„  +  tan  9j  -  tan  9^  ,      {2.k0) 

22.  O  d.  Vl  il 

1 

Ctn  9  ,  =  n^  -^^  (ctn  9  -  tan  9^^)  -  tan  9^^  .   (2.1fl) 

Now,  according  to  property  P  of  equation  (2.56),  no  generality  is 

lost  in  restricting  9„  to  the  range  -90°  <  9„  <  90°.   Assuming,  for 

n  11 

the  moment  that  this  restriction  has  been  made,  we  readily  verify, 
by  referring  to  Figures  (5)  and  (^) ,  that 


0-22 
^2 


and 

^24 


^2 


=  1 
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Introducing  the  acute  angle  of  reflection^ 


«  -  0^2  .   if  ©22  ^  °  ^ 


©22  =  (2.1^0) 

j^-(:t  +  922)  ,  if  ©22  <  0  , 


into  equation  (2.58)^  we  then  obtain  the  formulas 


Gtn  0^„  =  ctn  9^  +  2  tan  9^  ,  (2.4l) 

dd.  d  n 


Ctn  92i^  =  -^   ctn  Q^   +  {^  -   l)  tan  9^^       (2.1^2) 


which  are  valid  for  all  9^,  4  ±90. 

It  should  be  stressed  that  these  formulas  are  also  valid  for 
all  values  of  \Jf  +  90°j  i.e.,  ^  need  not  lie  in  the  plane  of  incidence. 
It  is  therefore  noteworthy  that  9  ^  and  9p,  depend  only  on  Q„,   the 
angle  that  H  ,  the  projection  of  H  on  the  plane  of  incidence,  makes 

with  the  normal.   This  Interesting  fact  seems  to  have  been  overlooked 

2 
by  V.C.A.  Ferraro  who  derived  another  somewhat  more  complicated  pair 

of  equations. 

To  relate  Ferraro 's  equations  to  ours,  we  refer  to  Figure  1  and 

note  the  following  relations  between  0,^  and  the  polar  angles  p  and  7 

n 

employed  by  him: 

H  cos  9^  =  H  sin  P 
P      H 

H  sin  9,,  =  H  cos  p  cos  7  . 
P      ^ 


f 
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These,  in  turn,  imply  the  relation 

tan  9„  =  ctn  3  cos  7 

n 
which,  combined  with  equations  (2.ifl)  and  (2. 14-2)  yields  Ferraro's 
results  -  apart  from  obvious  notational  differences. 

In  Figure  6(a)  we  have  plotted  9   against  0^  for  9^  =  0  , 

9  =  1+5°  and  9^  =~k^°-      These  are  typical  curves  in  the  parameter 
H  H 

classes  9„  =  0,  9_  >  0  and  9^  <  0,  respectively.   In  Figures  6(b)  and 
6(c)  we  show  plots  of  9  ,  versus  9  for  the  same  set  of  9  -values j 
however,  in  6(b)  t)  =  i)-  and  in  6(c)  ,  ti  =  l/k.   The  reader  may  wish 
to  check,  properties  P,  -  P  against  these  curves. 

The  angles 

9*2  =  ctn"^  [2  tan  9 J  (2.1^3) 

and 

1 

9*^  =  ctn"^  [(ti^  -  1)  tan  9 J  (2.1^) 

introduced  in  Figures  6(a)  -  (c)  are  the  maximum  positive  angles  of 
reflection  and  refraction  on  the  cu37ves-  of  the  class  9^-^  >  0,  t]  >  1  and 
the  mlniminn  negative  angles  of  reflection  and  refraction  on  curves  of  the 

class  9^  <  0,  Ti  >  1. 

n 

The  angles 

9^=   ctn"^[-2  tan  9j  ,  (2.1l5) 

and  .. 

9^"=  ctn"^[(Tl  ^  -  1)  tan  9  ]  ,  (2.46) 


26a 


^2=ctn-'[-2tan^J 


^^=ctn-'[-2tan^H] 


FIGURE   6(a) 
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02  =  ctn"'[(77^i)tane„] 


FIGURE  6       X 
(b) 


<^J=ctn-'[(75^-i)taneH] 


FIGURE  6 
(c) 
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are  the  positive  critical  angles  [See  footnote  Ik]    for  the  reflected 

and  refracted  waves^  respectively^  on  curves  of  the  parameter  class 

9  <  0,  T]  >  1  and  the  corresponding  negative  critical  angles  on 
n 

curves  of  the  class  0„  >  0,  ti  >  1.   Not  shown  in  Figures  6(a)  -  6(c) 

n 

are  the  angles  of  reflection  and  refraction  associated  with  angles  of 

incidence  9  which  produce  anomalous  reflection  and  refraction  as 

defined  In  the  previous  section.   These,  to  be  sure,  may  also  be 

computed  from  equation  (2.1t-l)  and  (2.ii-2).   Consider,  for  example, 

the  curves  of  Figure  6(a).   Here,  one  can  obtain  the  anomalous 

values  of  9   by  including  values  of  9  in  the  range  9   <  Qp  5  ^5 

J, 
when  9„  =  -^5°  and  by  including  values  of  9^  in  the  range  9^  >  9  >  -k^ 

when  9„  =  45  .   In  the  first  case  9^^  will  take  on  anomalous  values 
n  dd 

between  90  and  135  while  in  the  second  case  '^       will  take  on 
anomoloiis  values  between  -90  and  -155  •   In  both  cases  the  reflected 
wave  fronts  will  be  fo\md  in  the  same  quadrant  as  the  incident  wave 
fronts. 

5.  Relation  of  the  Graphical  Procedure  to  Huyghens  '  Construction 
In  order  to  explain  the  relation  between  the  two  procediires 
it  is  enough  to  consider  the  reflected  waves  alone*  the  refracted 
waves  may  be  treated  in  an  entirely  analogous  fashion. 

Huyghens'  construction  may  be  described  as  follows:  At  the  instant 
the  incident  wave  front  reaches  a  point,  0  say,  of  the  interf ace  ^C/ ^ 
0  becomes  an  induced  source  and  emits  a  "spherical"  wave.   The  same 
is  evidently  true  of  each  point  of  the  curve  of  intersection  of  the 
incident  wave  front  and  oO  .   To  construct  the  reflected  and  refracted 
wave  fronts  at  time  t'  >  t  say,  one  simply  constructs  the  envelopes 
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'k 


to  the  spheres  emitted  at  all  times  t  <  t'.   In  Figure  T,   0'  represents 

the  trace  point  at  time  t'  and  0  is  the  position  of  the  trace  point 

one  second  earlier.   For  the  sake  of  concreteness  it  is  assumed  that 

0'  moves  with  the  constant  velocity  V  =  Vx  and  that  V  >  0.   The  distance 

of  0'  from  0  is  therefore  V.   On  the  assumption  that  r.  *„  and  0^ 

11       n 

are. 5;  0  and  k^   ,   respectively^  the  "sphere"  through  0  can  he  shown 

22 
to  he  the  surface  of  revolution  obtained  by  revolving  the  f  ,  A     and 

s  -curves  of  Figure  7  about  the  H-axis.   The  "spherical"  wave  fronts 

thus  traced  out  are  known  collectively  as  Fresnel's  ray  surface  through 

0.   The  W. 's  in  Figure  7  represent  the  three  possible  types  of  incident 

wave  fronts  which  can  have  the  trace  velocity  V,  and  the  W. . 's  represent 

•^  ij 

the  associated  family  of  reflected  wave  fronts.   The  W.  .  's  and  the 
dashed  extensions  of  the  W. 's  are  tangent  to  the  f  ,   A  and  s  -curves 
since  by  constmction  these  wave  fronts  are  the  envelopes  of  the 
"spherical"  wave  fronts  emitted  earlier  and  up  to  one  second  later  than 
the  "sphere"  shown.   It  shoiild  be  stressed  that  the  "hemisphere" 
through  0  which  lies  in  the  region  N  >  0  of  Figure  7  -  i-e.^  in  ^  - 
indicates  merely  where  the  wave  fronts  emitted  from  0  would  have  reached 
were  there  no  discontinuity  across  *C7 .   The  "hemisphere"  appropriate 
to  IK     does,  in  fact,  differ  from  that  shown  when  Jj   is  a  surface  of 
discontinuity  and,  as  in  geometrical  optics,  can  be  employed  to  construct 
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arametric  equations  of  this  surface  are  given  in  both  the  paper  and 
report  of  Reference  1. 


T?he  Alfven  wave  fronts  evolving  from  0  consist  of  the  two  points 
labelled  A  in  the  figure.   For  the  purposes  of  Hioyghens  '  constr 

these  two  points  should  be  regarded  as  small  spheres. 
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the  refracted  wave  fronts  by  Huyghens  '  procedure. 

We  shall  now  explain  the  relationship  of  Huyghens '  construction 
to  the  graphical  procedure  of  the  preceding  sections.   To  convey  the 
essential  elements  of  this  relationship^  it  will  suffice  to  show  how 

to  obtain  Figure  ka   from  Figure  7-   We  assert  that  Figure  ii-amay  be 

23 
obtained  from  Figure  7  first  by  constructing  the  pedal  curves   to  the 

curves  of  Figure  7  and  then  by  inverting  these  pedal  ciirves  with 

respect  to  the  unit  circle.   To  prove  this  we  first  observe  that  the 

f  ,   A  and  s  of  Figure  7  are  the  curves  of  intersection  with  the 

2k 
plane  of  incidence  of  the  "sphere"  through  0.    On  the  other  hand^ 

the  f^  ,  A  and  s_,  curves  of  Figure  k   are  the  curves  of  intersection  of 
12      3 

the  plane  of  incidence  with  the  surface  of  wave  normals  (b  serving  as 
the  unit  of  speed) .   Since  the  surface  of  wave  normals  can  be  obtained 

from  the  pedal  surface  of  the  corresponding  Fresnel  ray  surface  by 

25 
inverting  with  respect  to  the  unit  sphere^   it  follows  that  the  curves 

of  Figure  i)-(a)  may  be  obtained  from  the  pedal  curves  of  those  of  Figure  7 

by  an  inversion  with  respect  to  the  unit  circle. 

To  complete  the  picture  we  consider  the  pencil  of  lines  through 


The  pedal  curve  of  a  curve "G, whose  equation  is  x  =  x(^ ) ,   is  constructed 

as  follows:  Let  T(|)  be  the  tangent  to~C,  •      From  the  origin  0^  draw 

the  perpendicular  to  T(^ )  and  let  y  =  y(l)  denote  the  vector  from  0  to 

the  foot  of  this  perpendicular.   The  curve  traced  out  by  y  as  |  varies 

over  its  domain  of  definition  is  the  pedal  wave  of  i^  .   '~*' 

2.k 

This  statement  is  true  only  when  H  is  in  the  plane  of  incidence,,  as  is 

the  case  here. 

25 

For  a  proof;  see  e.g.,  the  lecture  notes  by  E.  K.  Luneberg  "Propagation 

of  Electromagnetic  Waves";  New  York  University^  New  York  (19'^9)  • 
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0'.   This  family  of  lines  includes  among  others  the  lines  representing 
the  six  wave  fronts.   Let  0  be  the  angle  that  the  nonaal  to  a  line 
of  the  family  makes  with  the  vertical  N-axls  through  0'.   In  Figure  7; 
for  example,  9.,  is  the  angle  that  the  slow  wave  front  makes  with  the 
K-axis.   Then,  it  is  easy  to  verify  that  (l)  the  pedal  cxirve  of  each 
line  of  the  family  consists  of  a  single  point  with  the  polar 

coordinates  (s,0)  =  (V  sin  Q,   9)  and  (2)  the  locus  of  these  points 

V 
is  a  circle  of  diameter  —  whose  equation  is 

s  =  V  sin  9  . 

If  we  invert  this  circle  in  the  unit  circle  -  i.e.,  if  we  replace  s 

by  — ,  we  find  that 
P 


pb  sin  9  =  ^ 


which  is  precisely  the  equation  of  the  line  i  of  equation  (2.19) 


The  present  discussion  applies  with  obvious  changes  when 

the  parameters  9„,  V/b  and  r  are  varied.   Thus  varying  9  has  the 

effect  of  rotating  the  f  ,  A  ,    s  -curves  rigidly  about  0.   Varying 

—  has  the  effect  of  moving  0'  closer  or  farther  from  0  and  at  the 
b 

same  time  of  governing  the  question  of  critical  reflection.   For  example,  if 
0'  falls  within  the  f  -oval,  then  a  fast  wave  is  not  reflected  as 
there  are  no  tangents  to  the  f  -curve.   Moreover,  if  0'  coincides  with  0'' 
we  obtain  a  family  of  wave  fronts  disposed  in  a  way  alluded  to  in  subsection 
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B-5;  specifically^  an  incident  slow  wave  front  at  0 '  '  induces  slow 
and  Alfven  reflected  wave  fronts  wMch  lie  in  the  same  quadrant  as 
the  incident  wave  front  itself.  Finally^  varying  r  has  the  effect 
of  changing  only  the  detailed  features  of  the  cirrves  of  Figure  7  so 
that  the  above  discussion  applies  with  little  or  no  change  for  all 
r  >  0. 

Note  that  increasing  it^   from  zero  has  the  effect  of  turning 

n 

the  H-axis,  the  axis  of  revolution  for  the  Fresnel  ray  surface^ 
out  of  the  plane  of  the  page.  When  \lf„  >  0^  Huyghens  '  construction 
therefore  involves  finding  the  points  of  tangency  to  a  family  of 
similar  "spheres".   It  is  easily  seen  that  these  points  of  tangency 
do  not  lie  in  the  plane  of  incidence  and,  in  fact,  are  not  even  coplanar. 
Nevertheless,  Huyghens'  construction  may  he  made  graphically.   One 
need  only  work  with  the  projection  on  the  plane  of  incidence  of  the 
locus  of  points  on  Fresnel 's  ray  surface  at  which  the  normals^  are 
parallel  to  the  plane  of  incidence.   Parametric  equations  of  these  curves 
may  be  obtained  from  the  ray-equations  ( 3.^1-5) -(3- ^7)  of  Section  III 
by  substituting  the  expressions  for  n  and  H^ given  in  equations  (3-5) 
and  (3-6). 

In  graphical  work,  the  method  presented  here  is  more  accurate 
than  Huyghens '  method  in  that  a  point  of  intersection  is  more  easily 
located  than  a  point  of  tangency.   Moreover,  since  the  angles  of  reflection 
and  the  (reciprocal)  disturbance  speeds  are  determined  simultaneously, 
it  yields  more  information.   The  avoidance  of  the  cusped  figures 
should  also  be  noted.   Unquestionably,  however,  Huyghens  '  construction 
has  more  direct  physical  appeal. 


-  52- 

III.  REFLECTION  AT  M   INFINITELY  CONDUCTING  RIGID  WALL 
A.   Summary  of  the  Basic  Equations 

Let  tS  =  J[  i'^)    d.enote  any  propagating  planar  wave  front  and 
n  the  unit  normal  at  each  point  of  ^  (t)  pointing  along  the  direction 
of  propagation.   A  region  will  be  referred  to  as  being  in  front  of 
or  behind  A  (t)    according  as  n  or  -n  is  directed  into  that  region. 
Let  Q  be  the  undisturbed  value  of  some  physical  quantity  and  Q  and 
€L   denote  the  values  of  this  quantity  in  front  of  and  behind  ^(t) . 
Here  and  hereafter,  we  shall  take  the  state  of  the  medium  behind  the 
wave  front  as  the  state  in  front  plus  a  small  perturbation.   Letting 
&Q  represent  such  a  perturbation  in  the  quantity  Q  ,  we  therefore  have 
Q^  =  SQ  +  Qq  .  (3.1) 

In  a  given  wave  0^^  may  be  regarded  as  the  "  umperturbed"  value  of  the 
physical  quantity  it  represents.   It  should  be  stressed,  however,  that 
owing  to  the  presence  of  other  perturbations,  Qq  might  be  the  sum 
of  Q  and  a  small  perturbation.   Whether  or  not  Qq=  Q,  or  Q  plus 
a  small  perturbation,  it  can  be  shown,  on  neglecting  all  but  the 

highest  order  terms  that  5H,  5u,  and  5p  as  Alfven,  slow  and  fast  modes 

26 
are  given  by  the  following  formulas : 


Alfven:   5H  =  -|r:  Hx  n  , 


5u  =  +  — Hxn  ,  (3.2) 


1  pc  •«  K, 


5p  =  5S  =  0  ; 


2S 

K.  0.  Friedrichs,  "Nonlinear  wave  motion  in  magnetohydrodynamics", 

Los  Alamos  Report  No.  2105  (written  1954  -  distributed  1957)-   See  also 

later  version  of  this  work  by  K.  0.  Friedrichs  and  H.  Kranzer,  Report 

No.  MH-8,  AEC  Computing  and  Applied  Math.  Center,  Inst.  Math.  Sci. , 

New  York  University,  New  York  (1958). 


-  35 


5H  =  -  -^  nx  (nxH)  , 


fE^ 


6u  =   +- 


T=Vt"')--#.^-'^'-^-''   "■" 


y^    \c^      /  '^  ~  yp; 


e    ,  \'    '^ 


6p  == p  \  -5 1 


"j^ 


\c 


In  these  formulas,  H,  H  ,  i_i;  p  and  S  are,  in  the  stated  order,  the 
magnetic  intensity,  the  component  of  magnetic  intensity  in  the 
direction  of  propagation,  the  specific  magnetic  inductive  capacity, 
the  density  and  the  specific  entropy  of  the  flow  all  in  the  absence 
of  any  distuirbance  whatsoever.   The  c's  are  the  disturbance  speeds 
appropriate  to  the  mode  in  question.   They  are  the  positive  solutions 
of  the  equation: 

"  ''   ^  ^   =  0.  (3.1^) 


(pc^  -  ^^  ) 


k        ,2  „2.  2    2  „  2 

pc  -  (pa  +  jiH  jc  +  a  M-H 


Specifically,  in  equation  (3-2),  c  is  the  positive  root  of  the  first 
factor,  while  in  equation  (3-3) ^  c  is  the  smallest  or  largest  positive 
root  of  the  second  factor  according  as  the  mode  is  slow  or  fast. 
With  n  and  H  expressed  in  terms  of  the  angles  shown  in  Figures  1  and  2 


I.e. 


n  =  cos  0  N  +  sin  9  x  ,  (3'5) 


-o 


H     =     H  cos  iV^Xcos   0^^  N  +   sin  0,^  x   )+  H  sin  \lf„z„,(3.6) 
v«v  n  H  '~'  n  *-0  ^il~-0' 

we   obtain  for  the   slow,   Alfven  and  fast  disturbance   speeds,   the  explicit 
expressions   simmiarized  in  equations    (2.^4-)    -    (2.6).      The  disturbance 
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speeds  along  n  are  related  to  the  speed  of  propagation  U  >  0  of  >^(t) 
along  n  and  to  the  trace  velocity  V  as  follows: 

WW  ^v^ 

Vn  =  U  =  u-n  +  c  ;  (j.?) 


\AV  VW 


(see  Section  II-A) .   The  upper  (lower)  signs  in  equations  (3.2)  and 
(5.5)  correspond  to  the  upper  (lower)  signs  in  equation  (5>7)'   Note 
that  there  will  be  only  one  speed  of  propagation  along  n,  specifically^ 
the  one  corresponding  to  the  upper  sign  in  the  above  equations,  unless 
the  projection  of  the  fluid  velocity  u  on  n  exceeds  the  disturbance 
speed  G. 

The  absolute  value  of  the  quantity  e  appearing  in  (5-2)  and 
(5.3)  measures  the  strength  of  the  discontinuity  at  the  wave  front 
of  the  mode  in  question.   For  |e|  to  be  a  true  measure  of  this  strength 
it  is  necessary  to  nomialize  the  mode  of  propagation  in  a  suitable 
.anner.   This  is  the  function  of  the  constant  factor  "j/e^  whose  choice 
for  each  mode  is  dictated  by  the  following  considerations:   In 
reference  1,  a  mode  of  propagation  was  defined  to  be  a  six- vector  of 
the  fona 


mi 


R  = 


1 

J  (3.8) 


('H.  '  2  5H   [V  ,  2  gH  ,  5u  ,  ou  ,  5u  ,  ^ 
\P  y    \P  z'   n-'   y'       z'     p 


here,  the  x-axis  is  assumed  to  be  directed  along  n.   R  represents 
the  total  discontinuity  in  the  medium  at  a  given  point  of  the  wave 
front.   A  total  of  six  modes  of  propagation  are  obtained  by  substituting 
the  appropriate  components  from  equation  (3-2)  or  (3-3)  into  the 
above  expression.   Now  it  is  natural  to  define  the  strength  of  the 
discontinuity  to  be  proportional  to  (R  •  R)  ,  the  Euclidean  length  of  R. 
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For  reasons  which  will  soon  be  apparent  it  is  convenient  to  choose 
Vp/s  as  the  constant  of  proportionality.   With  (-^pR  •  R)^  as  the 
definition  of  mode- strength,  it  then  follows  that  E^  must  be  chosen 
to  be 


E 


1  _ 


§  R^  Ri 


£i^,aM,   ,(,„.,2^   (j.g, 


II  1  ? 

if  |e|  is  to  equal  (^pR  •  R  )  .   In  equation  (3.9)  R^,  &1I,  5pi  and 
ffl*  denote  the  values  of  R,  5u,  Sp  and  5  H  obtained  by  setting  €  =  1 
in  equations (5. 2)  and  (3.3).   If  one  actually  inserts  the  expressions 
for  5i/.  5p^  and  8H^  in  equations  (3.2)  or  (3.3)  into  the  last  member 
of  equation  (3-9)^  one  finds  that 


Alf v^n : 


1^  =  y^  b  I  sin  6  I  , 


(3.10) 


Slow 
or 

Fast 


T^  =  Vp  b  Isin  sl  ' 


2w  -  (l+r) 


(w^-  r) 


=  yp^i 


w 


1 
2 


(w  -  cos  5) (w  -  (l+r) 


(3.11) 


In  these  equations  5  is  the  angle  between  n^  and  H.   Employing 
equations  (3.5)  and  (5.6)  it  follows  that 


WW  vvv 


|sinS|      =    Tl-cosiircos    (9-9..) 

■H.  H 


(3.12) 


The   quantity  w  is  defined  by  the  equation 


w     = 


2  > 


(3.13) 
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it  is  understood  here  that  c  represents  slow  or  fast  disturbance  speed 
according  as  the  mode  is  slow  or  fast. 

Note  that  E^  has  the  dimensions  of  energy  density;  E^  may  be 
thought  of  as  the  "energy  density  of  the  mode"  R^.   Aq  Interesting 
sidelight  concerning  such  energy  densities  is  the  following.   For  modes 
R  of  arbitrary  strength  [^see  equation  (5.8)]  ,   one  can  show^  by  means 
of  equations  {3.k)    and  (5.2)  or  (5.5)  [See  Appendix  IlJ  that 


1^2     1  a^Sp^  ^  1  =,1x2  r-K  DA 

gPBu   =  2  — ^■^2'^^^   '  (5.14) 


so  that  the  energy  density  of  the  mode  is  divided  in  equal  parts  between 
the  kinetic  energy  density  of  the  mode  on  the  one  hand  and  the  sum 
of  the  internal  and  magnetic  energy  densities  on  the  other.   From 
this  equation  it  follows  readily  that 


E  =  p5u^  ,  (all  modes)        (5-15) 


|e  =  |p5u^  =  ^m^  (Alfven  mode)  .    (5-16) 


Here,  E  denotes  the  energy  in  a  mode  of  arbitrary  strength. 

B.   Reflection  at  an  Infinitely  Conducting  Rigid  Wall 
1.   The  Boundary  Conditions 

We  now  imagine  that  the  upper  half- space  uC    in  Figure  1 
is  a  rigid  infinite  conductor  and  that  the  unperturbed  state  of  the 
fluid  in  the  lower  half-space  W   is  constant  with  the  fluid  velocity 
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u  =  0.   The  (planar)  incident  wave  Is  assumed  to  propagate  into 
the  unperturbed  fluid  and  the  perturhations  in  the  regions  between 
the  surfaces  of  discontinuity  -  i.e.^  the  wave  fronts  and  the  interface 
J^  -   are  assumed  to  be  constant.   At  the  interface  (^,   behind  the 
slowest  reflected  wave,  we  have  the  following  boundary  conditions 


N  •  &u*  =  0  , 
BC:  N  •  5H*  =  0   ,  (3.1?) 

N  X  5E*  =  nN  X  (H  X  Bu*)   =   0   , 

■X-  * 

where  5u  and  5H  represent  the  total  pertiirbation  at  the  interface. 
The  first  condition  expresses  the  rigidity  of  the  wall;  the  second 
the  conse2rvation  of  magnetic  flux  and  the  third ,  the  continuity  of 
the  tangential  electric  field  vector  across  ^.      If ,   as  will  here  and 
hereafter  be  assumed,  H  is  in  the  plane  of  incidence,  then  the  boundary 
conditions  split  up  into  two  sets,  namely, 

N  •  6u^   =  0   , 

BC  :  N  •  6H*   =  0   ,  (3-l8) 

p  -»    ~p         ' 

N  X  (5u   XH)   =  0   , 

and 

BC^:  N^x(5u_L*X  H)   =   0  .  (3.I9) 

The  first  set  evidently  involves  only  6u   and  6H   ,  the  components  of 
Su_ and  6H  in  the  plane  of  incidence,  while  the  second  involves  Su,  , 
the  component  of  u  normal  to  the  plane  of  incidence.   Now,  as  an 
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inspection  of  equations  (3.2)  and  (3-3)  reveals,  ou  and  6H  are 
normal  to  the  plane  of  incidence  when  the  wave  is  an  Alfven  wave 
and  parallel  to  the  plane  of  incidence  when  the  wave  is  slow  or 
fast.   Since  the  boundary  conditions  for  the  normal  and  planar 
components  are  uncoupled  it  follows  that  Alfven  incident  waves  can 
produce  only  Alfven  reflected  waves  and  slow  or  fast  incident  waves 
can  produce  only  slow  and  fast  reflected  waves.   Assuming  for  the 
moment  that  the  incident  wave  is  slow  or  fast,  we  note  that  there 
are  two  reflected  waves  and  apparently  four  boundary  conditions 
in  equation  (3'l8).   It  will  now  be  shown  that  the  conditions  BC 
of  equation  (3.I8)  reduce  to 

BC  5u  =  0  ,         if  H  •  N  ^  0  ,      (3-20) 

and  to 

BC  N  •  5u  =  0  ,      if  H  •  N  =  0  .      (3-21) 

P  vw         >~p  '  ~»       l«. 

First,  we  observe  that  equations  (3. 20)  and  (5.21)  follow 
directly  from  the  first  and  third  boinidary  condition  of  (3.17) • 
To  complete  the  proof  we  must  show  that  N  •  5H  vanishes.   To  this  end, 
we  note  the  following  relations  which  hold  for  any  propagating  wave 
[see  reference  26J 

n^  •  5H  =  0  ,  (3-22) 

c  5H  +  n  X  (6uxH)  =  0  . 

Replacing  c  by  V   •    n    [see  equation   (3-7) J   note  u  =  o] ,  we  find  that 

n  x(5HxV  +  SuxH)      =0.  (3-23) 
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This   relation,   in  turn,   implies  that 

5HxV  +   5uxH     =     0   ,  (5.2i+) 

since  each  vector  within  the  parenthesis  is  normal  to  the  plane  of 
incidence  when,  as  we  have  assumed,  H  is  in  the  plane  of  incidence. 
One  such  relation  applies  to  the  incident  wave  and  to  each  of  the 
reflected  waves.   If  we  add  these  relations  we  find  that 


5H*  X  V  X  6u*  X  H  =  0  .  (5.25) 

But  5u   X  H  vanishes  by  hypothesis;  hence  sH  x  V  vanishes,  which 
implies  the  desired  result,  namely  N  •  5H*  =  0. 


2.   Strengths  of  the  Reflected  Waves 

We  shall  employ  the  symbol  0.  with  the  single  subscript 

i,  to  denote  the  value  of  a  physical  quantity  Q  in  front  of  an 

incident  wave  of  type  i  and  oQ.  to  denote  a  perturbation  of  Q.  in 

the  sense  of  Section  A.   Similarly,  Q.   and  gO. .  will  denote  the 

corresponding  quantities  in  a  reflected  wave  of  type  j  due  to  an 

incident  wave  of  type  i.   We  shall  further  allow  i  and  j  to  take  on 

the  "values"  s(for  slow),  f(for  fast)  and  A(for  Alfven)  in  place  of 

the  numerical  values  used  earlier.   Thus,  for  example,  e  denotes  the 

strength  of  a  fast  incident  wave  and  e  „  denotes  the  strength  of  a 

sf 

fast  reflected  wave  produced  by  a  slow  incident  wave. 

Case  1:  Alfven  wave  incident  -  Assiome  first  H  •  N  4  0-  "^^   boundary 
conditions  for  this  case  are  J_see  (3.19)_] 

Sii^A  +  5%^  =  0  >  (5.26) 


UO 


since ^  under  the  assiomed  conditions,  only  an  Alfven  wave  can  be 
reflected.   Employing  equations  (5.2),  (3.5)^  (3-6)  and  (3. 10),  we 
readily  conclude  that 


"AAl 


'A' 


(3.24) 


so  that  the  strength  of  the  reflected  wave  equals  the  strength  of 

the  incident  wave  even  when  the  acute  angle  of  reflection  [see   equation 

(2.3T)j  does  not  equal  to  the  angle  of  Incidence. 

When  H'N  vanishes,  the  incident  wave  causes  no  additional 

MM    «W  ^ 

distiirbance  within  the  medium,  since  no  boimdary  conditions  have  to 
he  met  [see   equation  (3.19)J'   At  the  interface  behind  the  wave, 
however,  we  find  a  surface  current  density  proportional  in  magnitude 
to  |&H.|. 

Case  2;  Slow  or  fast  wave  incident  -  The  following  formulas  are  the 
starting  point  of  o\ir  discussion: 


5u 


-c 


1^ 


5u  =  e  Su 

cos(9-9jt)cos9^ 

2 

w 


(3.28) 


cos9 


W  + 


cos(9-9„)sinG„ 

Jl      n 


sin9 


w 


aso 


E   =  1^  b  |sin(9-9^) 


2w  -  (1+r) 


2 
w  -  r 


=  Vp  b 


(w^-cos^(9-9jj)(2w^-  1-r) 


2 

w  -  r 


1 

2 


1 

2 


(3.29) 


w 


2/k2 
c  /b   . 
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These  formulas,  which  apply  to  both  slow  and  fast  waves,  follow 
directly  from  (3.5),  (5-5)^  (3-6)  and  (3-11)  on  setting  iir^  =  0. 

Let  us  first  consider  the  special  case  in  which  H  •  N  =  0 


W^     t^V 


-  i.e.,  in  which  0„  =  90°-   Here,  only  one  boundary  condition 

n 

need  be  met,  namely  N  •  5u  =0  [see   equation  (3.2l)]J.   Using 
equation  (3.29),  one  can  easily  show  that  this  boundary  condition 
can  be  satisfied  for  all  angles  of  incidence,  by  a  reflected  wave 
hairing  the  following  properties: 

(1)  It  is  of  the  same  type  as  the  incident  wave. 

(2)  It  emerges  at  an  (acute)  angle  of  reflection  which  is 
equal  to  the  angle  of  incidence.   Moreover,  the  relation 

e,.  =  e^  ,  (3.30) 

holds, where  i  is  s  or  f  according  as  the  incident  wave  is  slow  or 
fast. 

Suppose  that  0  4  j:90°  -  i.e.,  that  E   which  is  in  the  plane 
of  incidence  is  not  in  the  plane  of  the  interface.   In  this  case, 
we  have  two  boundary  conditions  from  the  vector  equation  (3-20)  and 
(at  least  for  angles  of  incidence  which  are  sufficiently  small  in 
absolute  value) ,  two  waves  for  meeting  these  conditions  -  one  slow 
and  one  fast.   Employing  equation  (3.29),  we  may  express  the  boundary 
condition,  behind  the  slow  reflected  wave,  as 

e.   6u.  +  e. „  6uf„  +  €.8u.  =  0  ,         (5.51) 
IS  "^is     if  -■if    1-1        ' 

where  i  is  f  or  s  according  as  the  incident  wave  is  fast  or  slow. 
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From  this  relation  It  follows  directly  that 


is  -"is    "^if       i  ~a    -if  ' 


€.^  6uf„  X  5uf   =  -e.  5uf-  X  &uf-  , 
if  ~if    --IS       1  -i    --is  ' 


which  in  turn  implies  that 

(Suf  X  5uf^  •  M) 
-1    -if   - 


(3.52) 


^is     = 

s 

^^^Is^'^'^tf 

•   M) 

"if  = 

h 

(SH|-  X  Su^f^ 

•    M) 

I'o..   ■»         ,  ,     0..  t_ 

.«N 

(3.55) 


(5uf-  X  5u?-„  •  M)  ' 
^  -IS   ~if   *-' 


Here,  M  is  the  unit  vector  normal  to  the  plane  of  incidence.   Since, 
for  fixed  r  and  9  ^  90  ;  numerical  analysis  shows  that  the  numerators 
and  denominators  in  (3- 55)  vanish  at  isolated  values  of  9.,  if  at  all, 
one  may  conclude  that  multiple  reflection  is  the  rule  and  single 
reflection  is  the  exception.   Observe  that  the  vanishing  of  numerators 
or  denominators  in  (3' 35)  means  that  the  velocity  vectors  involved 
in  these  expressions  are  parallel. 

Combining  equations  (5-29)  and  (3- 3)  one  obtains,  with  the  aid 
of  standard  trigonometric  identities,  the  following  explicit  formulas 


for  the  e .  .  '  s : 


is      1  f  '  sf  ' 

(3.5^) 

^if  =  ^^i^s'/^if 


it3 


Here,   D       is  given  by 


D 


sf 


V.   w,  _(Ef  Ef_)2 
IS  if^    IS  if 


•  v.Jv.    -l)cos(9.    -9„)sin(0.  „-9„) 
if^    IS      ^         ^    IS      H'         ^    if     H' 


+  w^3(l-wJf)cos(0.^-9jj)sin(O.^-9j^)  \     j      (5.35) 


the  corresponding  expressions  for  D   (d   )  may  be  obtained  from 

S      1 

this  equation  simply  by  eliminating  the  subscript  f(s)  in  all  terms 
where  it  appears  as  a  second  subscript.   The  values  of  v. ,   w. .  and 
9. .  appearing  in  the  D's  may  be  obtained  by  the  graphical  procedure 

-'-J 

described  in  the  earlier  sections,  and  E. ,  E^   etc.  may  be  calculated 

1   is 

on  the  basis  of  these  results.   In  this  fashion  one  can  obtain  a 
complete  quantitative  solution  of  the  problem.   Without  resorting  to 
numerical  or  graphical  methods,  however,  it  is  still  possible  to 
arrive  at  some  idea  of  the  content  of  formulas  (j.jJl)  and  (5.35)  l^y 
choosing  special  angles  of  incidence  or  special  orientations  of  the 
magnetic  field,  or  both.   We  shall  examine  several  such  special  cases 
now. 

Let  us  begin  by  assuming  that  Q.    ^  0   and  9^  =  0  -  i.e.  ,  that 
the  incidence  is  oblique  and  that  H  is  normal  to  the  plane  of  the 
interface.   In  this  case,  one  might  suppose  that  the  boimdary  conditions 
can  be  met  by  a  single  reflected  wave  emerging  at  an  acute  angle  of 
reflection  equal  to  the  angle  of  incidence.   Assuming  that  a  slow 
wave  is  incident,  this  would  imply  [according  to  equation  (3'53)j  that 
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e    „   vanishes  .   However^  on  substituting  ©   =n-0  or-rt+0 
SI  ss        s  s 

[according  as  9  is  positive  or  negative]  into  the  expression 
for  e  ,   one  finds  that  the  numerator  will  in  general  not  vanish} 
the  vectors  &u'  and  6u'   in  eqviation  (5.52)  are  evidently  not 
parallel.   Here  too  then,  double  refraction  is  the  rule. 

Assume  now  that  0.  =  0  and  0^^  4  0  -  i.e.,  that  the  incidence 
is  normal  but  the  magnetic  field  is  oblique.   Setting  0.  =  0 
0.^  =0.   =  l80  in  equation  (5'53)  and  observing  that  v..  =  w, 

if     IS  ^         ^      '  '='         11      1 

and  E. .  =  E.  we  find  that 
li    1 


6.   =  e.  5^  1  =  s.f 

IS      1  is 


^If  =     h   \f  ^  =  ^'^ 


(5.56) 


where  &. .  is  the  Kronecker  delta.   Thus,  when  H  is  oblique,  a  normally 
incident  wave  produces  a  single  reflected  wave  of  the  same  type  and 
strength  as  the  incident  wave. 

Suppose  now  that  0^^  =)=  0.   Let  the  angle  of  incidence  approach 

0  as  a  limit.   Assume  first  that  the  incident  wave  is  slow.   It  then 

H 

follows,  on  taking  limits  in  equations  (3.5)^  that  the  resulting  slow 
incident  wave  reduces  to  an  ordinary  sound  wave  when  r  =  a  /b  is  less 
than  unity.   The  same  is  true  of  a  fast  incident  wave  when  r  >  1. 
When  r  =  1,  both  the  limiting  slow  and  fast  waves  move  with  the  speed 
of  sound,  but  the  associated  6u  »s  and  SH.'s  have  non-vanishing 

components  tangent  to  the  wave  front.   However,  a  suitable  linear 
combination  of  these  waves  furnishes  a  so\ind  wave  in  this  case  also. 
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Thus  for  all  values  of  r^  one  can  assume  that  the  incident  wave  along 
H  is  a  so\md  vave.   The  point  we  wish  to  emphasize  is  the  following: 
A  sound  wave  incident  along  H,  can  produce  a  pair  of  geniilne  hydro- 
magnetic  waves.   It  is  necessary  of  course  to  assume  here  that  the 
angle  of  incidence  is  small  enough^  when  the  incident  wave  is  slow^ 
in  order  that  the  fast  angle  of  reflection  he  real.   To  obtain  the 
actual  strengths  of  the  reflected  waves ^  one  simply  lets  9  approach 
0^  in  the  numerators  of  the  right  member  of  equation  (5* 5^)* 

5.  Energy  Relations 

Let  u^  p  and  H  be  time-independent  and  constant  over  a 
given  region  of  space  -  q    ,    say  -  and  let  u ' ^  p  '  and  H '  denote 
continuous,,  continuously  differentiable  perturbations  of  these 
quantities  in  U  •   On  multiplying  the  linearized  hydromagnetic  momentum 
equation  by  u'-  and  combining  the  result  with  the  remaining  linearized 
equations  one  obtains,  as  in  gas  dynamics,  an  "energy"  equation, 
namely, 


||-  + V  •  F'  =  0  ,  (5.37) 


which  expresses  the  relationship  between  the  "perturbation  energy" 
per  unit  volimie 


2 
E'  =  I  p  (u'2  +^  p.2  +  ^jj,2^  ^  (3_38) 


and  the  "flux  of  energy" 
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F'  =  a%'u'  +  [i(Hxu')  XH'  +  E'u  .        (5.39) 


Corresponding  to  equation  (5' 37)^  we  have  the  jump  relation 

-USE  +  SF-n  =  0  ,  '      (3A0) 

associated  with  a  discontinuity  at  a  wavefront  moving  at  each  point 
with  the  speed  U  >  0  in  the  normal  direction  n.   5E  and  5F  here  denote, 
in  general,  the  difference  in  the  energies  and  fluxes,  respectively, 
in  front  of  and  behind  the  wavefront.   Since  we  always  regard  a 

disturbance  behind  a  wavefront  as  a  perturbation  of  the  state  ahead 

28 
we  conclude,  using  equations  (j.jS)  and  (5- 39)^  that 

5E  =  "I  p  (6u^  +  a^Sp^/p  +  laSH^)  ,  (3-^1) 


6F  =  a^SpSu  +  |i(HX5u)  X  5H  +  SEu  ,         (3-^2) 

[see  equation  (3.l)J. 

In  Appendix  II,  It  will  be  shown  that 

6F   =  5E  s  j  =  s,A,f  ,     (3A5) 


27 

_F'  in  the  present  theory  plays  a  role  analogous  to  that  of  the  Poynting 

vector  of  electromagnetic  theory.   The  second  term  in  equation  (3- 39)^ 

namely  |j.(Hxu')  XH'  =  ExH,  is  in  fact  the  Poynting  vector. 

The  quantities  a,p^H  and  u^ in  these  expressions  will,  in  general, 
differ  from  the  values  of  those  in  the  totally  undisturbed  flow  by 
first  order  perturbation  terms.   However,  it  is  clear  that  we  may 
ignore  the  contribution  of  these  first  order  perturbations  to  a2,p, 
H,  and  u,  since  they  contribute  third  order  terms  to  6E  and  5F, which 
are  themselves  of  the  second  order. 
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Here^ 

J        J 

and  s  is  the  ray  vector  associated  with  the  given  mode  of  propagation. 

Specifically^  in  slow  waves 

b^r(H  /H^) 
s   =  Jii±^.a± ^ (H-Hn)  ,       (5.^5) 


—s     »-  —  s-"  —   .  -^      '    n' 


c 
s 


in  Alfven  waves 


Sa  =  Ji±  (^^/P)^  sgnCH  )H  ,  (5.1^6) 


and  in  fast  waves 


_     l3^r(H    /h2) 


Sf     =     ii±c^    +      i (H-H^n)    .  (3.1+7) 


c^C 


In  these  equations 


,2 


(1+r)^   -   l+r(H   •   nf/E^ 


.1 
2 


(5.H5) 
and   sgn(x)    is  the   sign  of  x.      For  all  ray  vectors   s  we  find  that 

s-n     =     u+c     =     U.  (5-^9) 

v~       w.  n  — 

As  is  the  case  with  propagation  of  electromagnetic  waves  in  crystalline 
media^  hydromagnetic  energy  is  guided  along  rays  which  are  not  in 
general  directed  along  the  normal  to  the  wave  front  [see   reference  1, 
Section  B-5] . 

Turning  now  to  the  problem  of  reflection  at  a  rigid  infinite 
conductor^  we  expect  that^  at  the  interface  "behind  the  last  reflected 
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29 
wave,  the  total  normal  flux  will  satisfy  the  equation 

since  neither  electric  nor  mechanical  energy  may  pass  through  the 
surface.  In  this  equation  the  summation  is  extended  over  all  the 
fluxes  5F  .  associated  with  the  reflected  modes  of  type  j  excited 
by  an  incident  mode  of  type  i  with  flux  bF^.  Combining  equations 
(5.i|-5)  and  (3.50)  ^^   fi^  obvious  manner,  we  find  that 

^.SF.  .s,  .-N  +  SE.s.-N  =  0,  (3.51) 

or,  alternatively,  that 

S~.p6u^.s. .-N  +  p5u?s.-N  =  0.  (3.52) 

But  with  our  choice  of  the  normalization  factor,  we  have 

2 
5E  =  e 

in  any  given  mode.   It  follows  directly  from  equation  (5. 52)  that 

Y^     e%     -1+   ef  s  -N  =  0  ,  (3-55) 

the  summation  being  extended  over  all  reflected  waves  excited  by  an 
incident  wave  of  type  i  and  of  strength  I e .  I  - 


'^Formally,  equation  (5.50)  is  a  special  case  of  equation  (3.^)  with 
both  u  and  U  set  equal  to  zero.   The  analogous  relation  in  electro- 
magnetic theory  expresses  the  vanishing  of  the  normal  component  at 
the  Poynting  vector  at  an  infinitely  conducting  interface. 
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With,  the  aid  of  equations  (j-'^?)  -  (3.^7)  one  can  easily  verify  this 

relation  in  the  simpler  cases  where  only  one  ways  is  reflected.   For 

example^  in  the  case  of  Alfven  wave  reflection^  equations  {'^.hS)    and 
(3-55)  yield 

2      2^ 

since  s..  =  -s..   The  truth  of  this  equation  now  follows  from  the  fact 

^AA         "A 


that    |e..|    =    le.l     [see   equation   (3.27)] 


'AAi    '  A' 

IV.  REFLECTION  AND  REFRACTION  OF  N(M4ALLY  nrCIDEMT  SLOW  ATO  FAST  WAVES 
AT  A  CONTACT  DISCONTINUITY. 

Suppose  that  ^ In   Figure  1  is  a  contact-discontinuity  surface 

separating  two  stationary  polytropic  ideal  gases.   As  in  Section II-B; 

let  it  be  assumed  that  the  entropy  in^  is  the  same  as  that  in  "^  so 

that  only  the  density  is  discontinuous  acrossotT.   When  the  excitation 

is  normally  incident  the  reflected  and  transmitted  waves  propagate 

along  N  and  -E,   respectively.   Thls^  coupled  with  the  equality  of  r 

pop 
and  r'  [see  equation  (2.32)]  implies  that  the  ratio  w  =  c  /h  is  the 

same  for  all  waves  of  the  same  type^  reflected  or  transmitted  and  has^ 

in  fact,  the  value 

1 
w^  =  |(l+r)  -fi  [(l+r)2  _  l,r  cos^]^   .    (l+.l) 

In  this  equation,  the  upper  or  lower  sign  is  intended  according  as  the 
wave  is  slow  or  fast. 

It  is  assumed  below  that  H  is  in  the  plane  of  incidence  and  that 
L  =|:  0,   The  boundary  conditions  at  (^  are  then: 


50 


fc  B^V  N  =  fy   5u'\  N  ^  &u^,  (1^.2) 


(i^.i+) 


21  a^Sp  =  ^  a'^5p  .  (1^.5) 


The  summations  in  the  left  members  are  taken  over  the  incident  and  all 
reflected  waves  while  those  on  the  right  are, taken  over  all  transmitted 
waves.   These  conditions  are  simply  statements  of  the  continuity  of 
velocity,  magnetic  intensity   and  pressure  disturbances  across<<7. 
oL/^incidently,  is  assumed  to  move  with  the  velocity  5ury  [see   equation 

We  shall  show  that  for  normal  incidence  the  boundary  conditions 
may   be  satisfied  by  a  single  transmitted  and  reflected  wave  of  the 
same  type  as  the  incident  wave.   To  this  end,  let  e.,  e  and  e'  denote 
the  strengths  of  the  incident  wave,  the  reflected  wave  and  the 
transmitted  wave,  respectively,  it  being  understood  that  the  e's 
refer  to  fast  or  slow  waves  according  as  i  =  f  or  s.   Combining 
equations  (5.3),  (3-5),  (3-6)  and  (3.11),  and  recalling  that  (l)  0  =  l80° 
or  0  according  as  the  wave  in  question  is  a  reflected  wave  or  not 
(2)  w  is  the  same  in  all  waves  of  the  same  type,  we  find  that  equation 
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At  normal  incidence  the  5H's  are  all  parallel  to  the  wave  fronts 

and  hence  perpendicular  to  N. 
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(J+.2)  implies  e  -  e '  +  e^  =  0  , 

{k.3)   implies  e  -  e'  +  e   =  0  , 

(k.k)    implies  e"b  +  e 'b '  -  e  .b  =  0  , 

(lj-.5)  implies  eb  +  e'b'-e.b  =  0^   i=sorf. 

The  four  conditions  {k.2)-{k.^)   thus  actually  imply  only  two  conditions 
which  may  be  expressed  as 


e  -  e '  =  -e.   . 

1 


ih.6) 


eb  +  e'b'  =  e.b  .  (i^-.7) 


The  solution  of  this  system  is 


1 

e 
^i 

b-b' 
b+b' 

2 
n    +  1 

1 

t 
e 

^i 

2b 
b+b' 

_     2n^ 

1 

ih.d) 


n  =  p'/p  .  ik.9) 


The  positions  of  the  incident^  reflected  and  transmitted  wavefronts^ 
along  the  N  =  z-axis,  are  determined  at  each  instant  by 

z^  =  wbt  ,  t  <  0  , 

z   =  -wbt  ,  t  >  0  ^  (14-.  10) 

Z  '    =    Wb  't   ;,  t  >  0   . 
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Plots  of  these  ec[uations  are  shown  in  Figures  8(a)  and  9(a).   In 
Figure  8(a),  it  has  been  assumed  that  p  =  i|p  so  that  the  transmitted 
wave  front  travels  at  half  the  speed  at  the  incoming  and  reflected 
wave  fronts.   In  Figure  9(9-)  it  has  been  ass\mied  that  p  =  -j-  p  so 
that  the  transmitted  wave  front  moves  with  twice  the  speed  of  the 
incoming  and  reflected  wave  fronts.   Figure  8(b)  and  9{^)    shows  the 
waveforms  of  the  ratio  5Qrp/5Q.  where  50^^  represents  the  total 
disturbance  of  any  quantity  Q  -  e.g.,  the  pressure  or  N-component 
of  velocity  -  and  50.  refers  to  the  disturbance  of  Q  carried  by  the 
incident  wave.   Evidently,  6Q_,  in  the  neighborhood  of  the  interface 
c>Ois  in  absolute  value  greater  or  less  than  the  50.  according  as  p 
is  greater  than  or  less  than  p. 

Corresponding  to  equation  (5.55)  we  have  the  equation 

which  expresses  the  conservation  of  flux  through  J^.   But  in  the 
present  case  it  follows  readily  from  equations  (5.^5)  ^^^   (5-^7) 


b   >~<     b    .*v      b 


wtj 


so  that  equation  (4.11)  reduces  to 

b'e'^  +  b£^   =   be^^   .  (4.15) 

It  is  easy  to  verify  directly  with  the  aid  of  equations  (4.8)  and 
(4.9)  that  this  relation  is  indeed  satisfied. 
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z  =  w  bt 
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FIGURE  8 


(a) 
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(b) 


FIGURE  9 
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V.   LIKEARIZED  TWO- DIMEI^ SIGNAL  STEADY  FLOW 
A.   Aligned  Fields"^  Case 

We  present  here  an  account  of  several  standard  steady-flow 
problems  for  the  purpose  of  illustrating  another  form  of  application 
of  the  theory.   We  consider  first  the  motion  of  a  hydromagnetlc  fluid 
which  is  deflected  through  a  small  angle  owing  to  an  abrupt  change 
in  inclination  at  a  boundary  wall.   The  basic  setup  is  shown  in  Figure 
10  (a).   The  unperturbed  fields  u  and  H  are  assiimed  to  be  parallel 
to  each  other  and  to  the  wall  [in  the  region  (t     of  the  figure^  and 
for  definiteness  to  point  in  the  negative  x-direction.   Thus  we  may 
write 


u   =  -u  X   ,  u  >  0  ^         (5.1) 

v-o        o~o  '  o  ' 


H   =  -H  X   .  (5.2) 

>*o       O-'O 


The  coordinate  axes  are  fixed  at  the  comer  0.   All  angles  are 
reckoned  from  the  positive  direction  of  the  vertical  a^is  -  the 
N-axis.   The  vector 

A  =  sin  ox  +  cos  a  N  (5- 3) 

represents  the  unit  normal  on  the  inclined  portion  of  the  wall  h.- 
It  will  be  supposed  that  the  deflection  of  the  fluid  is 


■^  The  reader  is  referred  to  the  paper  "Compressibility  Effects  in 
Magietoaerodynamic  Flows  Past  Thin  Bodies"  by;  J.  E.  McCune  and  E. 
L.  Resler,  Jr.  for  a  comprehensive  account  of  this  subject  from  a 
different  point  of  view. 
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/ 
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(c) 


(d) 


(e) 


FIGURE   10 
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effected  by  one  or  more  weak  hydromagnetic  disturbances  of  the  type 
discussed  earlier  [see  SectionlH-A] .   Since  the  wavefront  normals 
n  and  all  disturbances  6u,  &H  etc.  of  the  state  in  front  of  the 
wave  are,  by  symmetry,  parallel  to  the  plane  detemilned  by  H  and  N, 
it  follows  that  Ali"ven  waves  may  be  eliminated  from  consideration 
at  the  very  outset j  for  in  Alfven  waves  5u  and  5H  are  both  normal 
to  the  plane  of  n  and  H  and  hence  to  the  plane  of  N  and  H. 

Employing  equation  (2.3),  we  find,  since  U.  =  0,  that 

u  X   •  n .   =  u  sin  9   =c,    O<0<n, 
o-o   *~i      o     i      i  '      —  i  - 

i  =  s,f,       i^.k) 


or  alternatively  that 


p.b  sin  ©i  =  jT  ^  ^5.5) 


"b 


where 


Pi  =  ~   ^ 


(5.6) 


and 


\ 


u/h^   ,    b^  =  l^jTV^  .    (5-7) 


o'  o       o    '   o 


In  these  equations,  n.  is  the  nonnal  to  the  i-th  wave  front  which 
points  upstream  into  the  region!^  ,  c.  is  the  disturbance  speed  in 
the  direction  n. ,  and  0.  is  the  angle  that  u.  makes  with  N.   Note 


that  9.  is  permitted  to  exceed  90  . 


-  55 


As  In  SectlonH,  to  determine  the  0.  's  and  c.  's  we  must 

1        1 

find  the  intersection  of  the  line 

£(M^):  pb^sin  ©  =  ir  (5.8) 

with  the  waves  obtained  by  plotting  p  and  p.  against  the  angles 
9  and  0  respectively  on  the  same  polar  diagram.   In  Figure  10(b) 
we  have  plotted  only  the  upper  parts  of  these  curves j  the  lower 
parts  may  be  obtained  by  reflection  of  the  upper  parts  in  the  H  -axis. 

Let  it  first  be  assumed  that  the  undisturbed  particle  velocity 
exceeds  or  is  equal  to  the  fast  disturbance  speed  in  the  H  -direction 


i.e.;  that 


f  -i     1 

M~     <  man  'r        }   ^  '{       >  r=a/b  .      (5-9) 


Then,  as  a  mere  glance  at  Figure  10(b)  reveals,  •?(M,  )  Intersects  the 
slow  branch  only  at  p  =  oo .   This  implies  that  9=0,  which  means 
that  the  slow  shock  front  is  parallel  to  the  flow.   Evidently,  of 
the  two  possible  shocks,  only  the  fast  shock  is  available  for 
effecting  the  deflection. 

Similarly,  when  M,  satisfies  the  relation 

,  r  ^,   l\     <     M^-^  <  [(l+r)/rj^  ,        (5-10) 


max 


only  the  slow  branch  is  intersected  by  i(M,  )  and  hence  only  a  slow 
shock  is  available  for  effecting  the  deflection.   When  K     satisfies 
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the  relation 


or 


1 
min  .1  r  ^ ,  1  [   <   K, " 


M^"^  >  [(l+r)/r] 


<  max 


(5.11) 


(5-12) 


neither  the  fast  nor  the  slow  shock  is  available  to  effect  the 
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deflection J   another  type  of  solution  must  be  employed. 

In  the  discussion  which  follows ^ the  subscripts  s  and  f  will 

be  omitted  with  the  understanding  that  any  given  formula  refers  to 

a  fast  wave  solution  when  NL  satisfies  the  relation  (5-9)  and  to 

a  slow  wave  solution  when  M,  satisfies  the  relation  (5. 10). 

The  remaining  part  of  the  discussion  is  based  upon  the  formulas 


cos  9  N  +  sin  9  x 


5H 


6u. 


5p 


eH 


^f^ 


ec 


V? 


2 
=  a^  5p 


°  '  sin  9  cos  9  W  -  (l  +  sin  9)  x  I   , 


w^-l 

cos  9  N  +  T—  sin  9  x 

2       ■«-'o 


w 


e    2 


y? 


a  p 
00 


1  - 


2  -, 
sin  9 


w 


(5.15) 
(5.1i^) 

(5.15) 

(5.16) 


when  II  satisfies  the  relation  (5.9)  or  (5.IO)  the  flow  is  hyperbolic; 

when  M,  satisfies  either  (5. 11)  or  (5-12)  the  flow  is  elliptic.   These 

hyperbolic  and  elliptic  flow  regions  are  precisely  the  hyperbolic  and 
elliptic  flow  regions  described  by  J.  E.  McCune  and  E.  L.  Resler,^  as 
can  be  shown  by  making  obvious  notational  changes. 

See  reference  h. 


57 


1+  r  +  [(l  +  r)' 


kr   sm  9J 


r  = 


2/T.2 
o'    o   ' 


(5.17) 


w 


(5.18) 


E   =  Vp  b   I  cos  gI 


gy  -  (l  +  r) 


2 

w  -  r 


1 
2 


(5.19; 


and  the  boundary  condition  at  the  (rigid)  wall 


u^  •  A  =  0 


(5.20) 


The  first  set  of  fonnulas  is  obtained  from  the  corresponding  set  in 
Section n-A  by  setting  0  =  90°  in  the  formulas  of  that  set.   Q^ 
here  refers  to  the  value  of  a  quantity  Q  on  the  wall  side  of  the 
shock  1$   [region"^   of  the  figure]  and  Q  to  the  value  in  the 
unperturbed  region"^  .   In  equation  (5. 17),  the  upper  sign  applies 
to  the  fast  wave  and  the  lower  to  the  slow  wave. 

From  the  shock  relations  for  steady  hydromagnetic  flow  it 

follows  that  if  H   is  collinear  with  u  ,  then  H,  is  collinear  with 
-^o  *~o      -^-l 

u  so  that  we  may  replace  the  boundary  condition  (5-20)  by  the 
equivalent  condition 


H. 


A 


(H  +  5H)  •  A  =   0 


(5.21) 


Note  that  it  is  immaterial  whether  the  wall  is  infinitely  conducting 
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or  not;  for  the  vanishing  of  both  u  •  A  and  H  'A  is  equivalent  to 
the  vanishing  of  the  tangential  electric  field  AX|a(H  xu  )  and  of 
u  •  A^  at  the  wall  CL. 

Combining  equations  (5-2)^  (5-3)  (5'13)  and  (5-l^)  we  find, 
retaining  only  first  order  terms ,   that 

£H 

-H  a  + sin  9  cos  9=0.  (5.22) 

33 
This  equation  and  the  equations 


sin  9  =  —  =  -   ;     ,  (5-23) 

o 


iv^ 


cos  9  =  +    '^    -^      ,  (5.2^) 


M 


lead  directly  to 


.V?-2 


+a  r  E  M 

e  =  ——-——  (5.25) 


^f^r-^ 


for  the  value  of  the  weak  shock  strength.   The  value  of  the  pressure 
behind  the  shock  is  therefore  [see  equation  (5.l6)J 

P  u  (a/c)  a   / 


The  minus  sign  in  equation  (5.2^4-)  is  associated  with  values  of  9 


exceeding  ^Q  . 
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similar  formtilas  may  be  obtained  for  5H  and  5u.   When  H  approaches 
zero,  both  b  and  mT  =  b  /u  approach  zero  and  c  approaches  a  ,   the 
speed  of  soimd  in  the  region  of  undisturbed  flow.   As  is  to  be  expected, 
if  the  solution  associated  with  the  lower  (minus)  sign  is  ignored, 
the  resulting  expression  reduces  in  this  limit  to  the  well-known 
gas  dynamical  expression. 

The  solution  associated  with  the  minus  sign  in  equation  (5.26) 

merits  special  attention.   Let  us  first  assume  that  the  deflection 

of  the  fluid  has  been  effected  by  shocks  -  weak  shocks,  to  be  sure, 

but  shocks  nevertheless.   Let  us  first  assume  that  lyL  >  1.   In  this 

case  the  shock  must  be  fast  [cf.  (5-9)]  an<i  we  must  choose  the  "plus" 

sign  in  equation  (5.26)  to  insure  that  the  shock  is  compressive. 

This  implies  that  the  shock  front  is  tilted  back  toward  the  wall. 

[See  Figure  10(a)J  .   Now  suppose  that  R  <  1.   Then  in  virtue  of 

equation  (5. 10)  the  shock  is  slow  and  because  the  factor  1  -  M" 

in  equation  (5-26)  is  negative,  we  must  now  employ  the  expression 

for  p^  with  the  minus  sign  in  order  to  insure  that  p,  exceeds  p  . 
1  °  1         o 

The  choice  of  the  negative  sign  implies  that  the  slow  shock  front  is 
tilted  upstream.-^ 

This,  at  first  startling  fact,  is  best  understood  by  considering 
what  happens  when  a  is  large,  and  strong  shocks  are  required  to 
effect  the  change  of  direction.   In  this  case  too  it  is  a  fact  that 
the  magnetic  field  must  be  horizontal  in  front  of  the  shock  and  parallel 
to  the  wall  behind.   From  this  and  the  continuity  of  the  normal 


Already  observed  by  J.  E.  McCune  and  E.  L.  Resler,  Jr. 
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component  of  H  across  the  fronts  it  follows  immediately  from  obvious 
geometric  considerations  that  (l)  the  tangential  component  H   on 
the  wall  side  of  a  shock  must  exceed  the  tangential  component  H 
on  the  upstream  side  if  the  shock  is  tilted  back  toward  the  wall 
and  (2)  the  relation  H   <  H   if  the  shock  is  tilted  upstream. 
As  is  well  known,  the  first  condition  is  satisfied  by  fast  shocks 
only,  and  the  second  by  slow  shocks  only. 

It  has  thus  far  been  assumed  that  a  is  positive.   If  a  is 
negative,  it  appears  from  ccnsiderations  similar  to  those  above  that 
the  deflection  of  the  stream  can  be  effected  by  a  slow  shock  tilted 
toward  the  wall -vAienll  <  1  and  by  a  fast  shock  tilted  upstream  when 
R^  >  1! 

If  we  remove  the  restriction  that  the  deflection  of  the  stream 
be  effected  by  shocks  only  and  regard  the  formulas  (5.1^)  -  (5 -16) 
as  giving  the  variation  of  the  magnetic  intensity,  the  velocity  and 
pressure  across  a  weak  simple  wavej  it  then  appears,  since  such  waves 
may  be  expansive,  that  the  deflection  can  be  effected  by  a  slow  expansive 
simple  wave  when  M,  <  1  and  a  >  0  and  a  fast  expansive  simple  wave 
when  M,  >  1  and  a  <  0.   In  this  way  it  is  possible  to  get  rid  of 
shocks  facing  upstream  if  such  shocks  are  found  to  violate  other 
requirements  on  the  flow. 

The  expression  for  p  in  equation  (5.26)  and  the  corresponding 
expressions  for  5H  and  5u  apply  with  obvious  changes  to  nose-on  flow 
past  a  symmetrical  wedge  with  wedge  angle  2a  [See  Figure  10(c)] .  Indeed, 
the  expression  for  p  in  (5.26)  already  gives  the  pressure  above  and  below 


-  61  - 

the  wedge  J  no  modifications  are  necessary.   Moreover  by  applying 
formula  (5.26)  at  each  comer  of  the  symmetrical  wing 'It/',  made  up 
of  linear  segments  [see  Figure  10(d)] ,   we  can  obtain  the  pressure 
distribution  at  each  point  of  the  wing.   In  fact,  if  we  imagine  the 
number  of  segments  to  increase  without  number  and  to  approach  the 
continuous  profile  shown  in  Figure  10(e)  then  the  expression  for 
p  in  equation  (5.26)  gives  the  pressure  on  top  and  on  bottom  of 
the  wing  at  a  distance  x  from  the  leading  edge.   Here,  a  must  be 
regarded  as  a  function  of  x. 

The  above  analysis  is  also  applicable  to  flow  past  a  thin 
plate  having  an  angle  of  attack  a.  In  this  case  the  pressure  of 
each  point  of  the  plate  is 


Pt   =  P  + 
1      o 


o 

2p  u  (a/c)a 
00' 


V^TT 


^^J 


(5.27) 


The  additional  factor  of  two  Is  due  to  the  fact  that  the 

pressure  forces  below  the  wing  add  to  that  above.   The  lift  coefficient 

due  to  pressure  force  alone,  namely, 


2  o  o 


is  therefore 


C^(M)   =   ^(^/"^^    f  1  -  -^  ^  (5.29) 

fV?  -  1 
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B.   The  Crossed-Fields  Case 

The  crossed  fields  case  can  be  treated  by  methods  analogous 

to  those  of  Section  A^  provided  M.  is  sufficiently  large.   Consider^ 

for  e.xample^  the  case  where  the  magnetic  field  is  parallel  to  N  and 

hence  perpendicular  to  the  direction  of  flow  -  u  x  .   In  this  case 

o-o 

the  wave  front  angles  are  obtained  from  a  figure  like  that  of  Figure 
11.   Obseirve  that  according  as 


\' 


> 


fr  +1 


or 


\ 


-1    .     1 


< 


1/r  +1 


one  or  two  disturbance  waves  are  available  for  meeting  boundary 

conditions  at  the  wall.   The  first  relation  defines  the  so-called 

35 
"hyper-liptlc"  region,   the  second  relation,  the  hyperbolic  region. 

The  method  of  Section  A  applies  in  general  only  to  the  hyperbolic 

region. 

In  the  above,  in  keeping  with  the  assumption  that  the  problem 
is  a  two-dimensional  one,  it  was  implicitly  ass\imed  that  -u  x  ,  N 
and  the  prevailing  magnetic  field  H  were  co-planar.   When  H  is 
not  in  the  plane  of  -u  x  and  |[,  then,  as  the  reader  may  readily 
verify  for  himself,  there  will  be  a  hyperbolic  region  of  flow  in 

which  Aljfven  as  well  as  slow  and  fast  waves  are  available  and  in  fact 

necessary  for  meeting  the  boundary  conditions. 
_ 

See  reference  51 


-  62a  - 


r  = 

•5.  eH=o« 

i 

>  N-  axis 

^H  =  0° 

1 



1 

— 

s 

^— 

1 

A 

/^ 

^ 

II  ,  ,                        ....      II 

—   Hyper-liptic  -► 

( 

A 

^-"Hyper- 

. .    II 
iptic    - 

► 

Vr+i 

/ 

fV 

^ 

/ 

A 

1 

s — 



. 

4— 

-K 

Hyp 

erbc 

)lic 

FIGURE  II 


-  63  - 

APPENDIX  I 


We  prove  that 


V  •  n.   =  U.    ,  •     (1) 


V  ■  n.  .  =  U.  .   ,  (2) 


and  that  n. ,  n. .  and  W  are  coplanar.   Our  proof  does  not  require 
the  wave  fronts  to  "be  planar. 

We  "begin  by  expressing  the  equations  of  all  wave  fronts  in 
the  form 


W(x)  -  (t-t^)  =  0  .  '  (5) 


Defining  p  by  the  equation 

p  =  vw/|vw|  ,  (h) 

we  see  that 

n  =  p/p   .  (5) 

It  is  also  easily  proved  that  the  speed  U  of  the  wave  front  along 
n^  is  related  to  p  by 

U  =  p-^  .  (6) 

At  each  instant  the  incident  and  scattered  wave  fronts 
intersect  the  interface  (lO  in  a  common  curve  W(t)  -  the  so-called 
trace  -  which  sweeps  but  a  portion  of  i/J  .      Let 


-  Gk  - 

be  the  equation  of  a  trace  for  each  t.   Without  loss  of  generality, 
it  may  be  assumed  that  for  fixed  i,   the  curve  y  =  y(|,t)  is  orthogonal 
to  the  family  of  traces  ol  (t) ,  t  >  t  ,    so  that  we  may  write 


dt 


V  =  ^        .       .  (7) 


I'low;  since  ^(t)  is  on  all  wave  fronts,  it  follows  that 

Wi[y(i.t)]    =   w..[^(^,t)]    =    t-t^  ,  (8) 

for  each  i  and  j.   Differentiating  first  with  respect  to  |  and  then 
with  respect  to  t,  we  find  that 


£i'  ^   =   £ij*  1  =    ^    '  (10) 

which,  using  equations  (5)  and  (6),  become 

n.  •  -rf  =  n.  .  x^  =   0   ,  (9) 
n.  -V     n.  .-V 

^  =  T^  =  1  ■  (10) 

Equations  (l)  and  (2)  now  follow  immediately  from  equation  (lO). 

Moreover,  as  IT  is  proportional  to  the  cross-product  of  ^  and 


Ti 


V  =  -^  it  follows  from  (9)  and  (lO)  that 


N  X  n.   =  N  X  n.  .   ,  (ll) 

which  implies  that  the  n. .  are  in  the  plane  of  W  and  n.. 
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APPENDIX  II 


We  sketch  here  a  derivation  of  the  relations 


1  _  2     1  2  5p"   1  _„2 
2^     =  2    p~  +  2^5H 


1  _  2   1  a  „  2   1  ^„2 
5E  s  2  p&u  +  2  ~  ^P  ^   2*^ 


p5u 


5F  =  a  5p5u  +  Li(Hx5u)x  5H  +  5Eu  =  5Es  , 


(1) 

(2) 
(3) 


which  hold  for  all  modes  of  propagation.   In  (5)  ^  denotes  the  ray 
vector  at  each  point  of  the  vave  front  of  the  mode  of  propagation 
in  question.   We  shall  also  show  how  to  derive  the  formulas 

2 


}e^   =     h  Vp  |sin  6| 


2w  -  (1+r) 


2 
w  -  r 


w 


by^  r 


w 


(w^- cos^5)(2w^-  (l+r)) 


T  1 
2 


2    2/2 
w  =  c  /a  , 


for  slow  and  fast  modes  and 


y?=TP 


b   sin  5 


(5) 


for  Alfven  waves. 

It  will  be  enough  to  give  the  derivation  for  the  fast  and 
slow  waves  since  the  procedure  for  Alfven  waves  is  quite  analogous 
and  in  fact  much  simpler. 
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We  begin  by  observing  that 


f   2        2w2,2v  h       2,    2^2,  2^2 

(c     -aj(c    -b    )   =   c    -c    (a+bj   +  a  d 
n  n  n 


2^>.2     -2. 
=   c   (b    -  b      ) 
n 


(6) 


c  lap"    (H-H^n)         , 


in  virtue   of  the  equation 


h        ,2.       ^2s    2  2^2 

i      -   ( a     +  b   )  c     +  aHD       =   0 
^  n 


(7) 


|_cf.    equation   (3.i|)J-      Employing  equation   (3-3)    anti   (6),  we  find 


5u 


2  2 

€    C      < 


2   2 

€     G 


b2\2        b   ^      (H-Hn)2 


n* 


/ 


H 


h 


2\  2 


■^ 


V 


n 


b 


1.  2    2 
b     c 
n 


1- 


---2 


(8) 


2      /       b 

2  2 

c        \        c 


c   -  a  b 


Performing  similar  calculations,  we  can  readily  relate  all  quantities 

2 
appearing  in  the  definitions  of  5E  and  SF  to  p5u  .   The  final  results 

are: 


p     c  (c  -  a)     2 
c  -a  b 


(9) 
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2^2       a^Cc^-b^) 
a  5p  n  _   2 

-f-  =      1,       P     2        P^^        ' 
c    -  a  D 


n 


(10) 


c5ff,2  _      2^ 
n(HX5u)X5H     =     +  -ff — „     _  %5u  n      , 


La^ 


n 
2„  „         a^c    L    2   ^  2  \   „  2 


a  5p5u  =  — ^-2-  Tn  h"  '  ^  S.   P^^ 


^-a^ 


(11) 


(12) 


n 


Adding  the  first  twiD  of  these  equations,  we  obtain  equation 
(l)  from  which  equation  (2)  follows  immediately. 

To  prove  (3),  we  observe  that  the  ray  equations  (3A6)  and 
(3.47)  may  be  expressed  as  follows: 


a^  ^  (H  -  H  n) 


u  +  cn  + 


n      n** 


H  (+  c  C  ) 
n^ 


(13) 


Here,  (l)  the  upper  and  lower  sign  in  the  denominator  of  the  last 
term  applies  when  the  rays  are  "slow"  rays  and  "fast"  rays, 


respectively,  (2)  C  is  defined  by 


C   =  I 


2    2       2.  2 
(a  +  b  )   -  ka-b 
n 


1 
2 


(Ih) 


[cf.  equation  {^.kQ)']    and  (3)  c  is  related  to  C  by 


2      2    2—2 
2c   =  a  +  b  +  C 


(15) 


as  may  be  seen  by  solving  equation  (7)  for  c  .  With  the  aid  of 
equation  ( 7) ^  we  find  that 
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o2  2 
+  C  c 


a  2      h 

a  b   -  c 


n 


(16) 


from  which  it  follows  that 


s  =  u  +  en  + 


a^  ^g(H-  H  n) 
n       n"' 

H  (a^  2_  ^1.^ 
u    n 


u  + 


^  ^   _2^2 


c  -  a 


A^^ 


,     aT3  H 
4       n  '^ 


(17) 


Now  adding  equations  (ll)  and  (12) ,   we  find  that 


a  5pSu     +h(HX6u)x5H     =     +  -^ r — ;r 

c  -  a  D 


a^2  - 

+       n  „ 

en-  — H 

n 


P6u   .   (18) 


Comparing  the  right  memher  of  this  expression  with  the  last  term  in 
equation  (17) ;  we  conclude  that 


a  5p5u  +  u(Hx6u)x6H  =■  p5u  (s,-u)   ;, 


(19) 


which^  using  equation  (2)^  leads  directly  to  equation  (j). 


To  prove  (h)   we  must  first  calculate  Rn'R-,'  The  simplest 
procedure  is  to  imagine  that  at  a  given  point  of  the  wave  front  a 
special  coordinate  system  has  been  introduced  with  the  x-axis 
directed  along  n  and  the  z-axis  taken  along  the  direction  of  SH^which 
is  tangent  to  the  wave  front.   It  is  then  easy  to  show  that 


K^.K^  =  .^..„^5  [flp^]      M^  .  A'-^^^        .    (20) 


T 
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where  5  is  the  angle  between  H  and  n.   Now  from  {&) ,   it  follows  that 


.^2  .  2. 


/  2  ,  2v     CD  sin  5 

^^  -\   ^    =    2   2 
c  -  a 


(21) 


2    2 
Substituting  for  (c  -  b   )  in  (20)^  we  then  find^  after  simplifying, 

n 


that 


\-«l=  2^%^ 


r  r     h 


(    ^       2,  2 
(c  -  a  D 
n_ 

2   2 
c  -  a 


(22) 


which,  employing  (7);  becomes 


R«R   =  2b  sin  5 


2c^-  (a^  +  b^) 

2    2 
c  -  a 


(25) 


Alternatively,  employing  (6),  R-. 'K-.  '^^^7   ^s  expressed  as 


\-\  -   2 


[2c2-(a2  +  b2)-|  [c^-b/j 


(2l|) 


Equation  {h)    now  follows  from  Equation  (25)  and  (24)  on  multiplying 
p/2,  introducing  w  for  c  /a  and  finally  extracting  roots. 
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